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SUMMARY 

The  McDonnell  Douglas  fatigue  crack  growth  model  CNTKM8  has  been  used  to 
predict  crack  growth  in  CCT-specimens  with  through  cracks  under  four  different  load 
sequences,  including  a  flight-by-flight  sequence,  and  two  stress  levels.  The  predictions  were 
compared  against  experimental  data.  The  results  ranged  from  fairly  good  for  the  more- 
structured  sequences,  to  rather  poor  for  the  more-random  sequences.  Furthermore,  the 
experimental  trend  of  the  structured  sequences  giving  longer  crack  growth  lives  than  the 
random  sequences  was  not  predicted  correctly.  An  enhanced  model  was  proposed,  which 
was  also  tested.  The  predictions  of  this  model  were  more  consistent  with  the  experimental 
data,  and  the  above-mentioned  experimental  trend  was  also  predicted  correctly. 

This  report  contains  a  detailed  description  of  both  models,  including  the  equations 
that  are  relevant  for  the  case  of  a  CCT -specimen  with  a  through  crack,  and  of  the  tests  that 
they  were  subjected  to. 
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NOTATION 


The  following  list  of  symbols  is  used  throughout  this  report  Where  possible  the  usual  notation 
as  found  in  the  international  literature  has  been  sustained.  Since  most  of  the  equations  have  been 
obtained  from  the  listing  of  the  computer  program  that  incorporates  the  Me  Air  model  [6],  however,  it 
sometimes  was  mote  convenient  to  use  that  notation. 


a 

aoi 

check 

da/dN 

cffkr 

fr 

Kc 


Kmb, 


^^in.ol 


Nd 

R 


Wol 

*fr 


half  crack  length  in  centre-cracked  tension  specimen 
half  crack  length  at  time  of  overload 

parameter  that  is  used  in  test  to  decide  whether  should  be  updated 
crack  growth  rate 

proportion  of  opening  stress  to  overload  in  overload  cycle  for  R  =  0 
ratio  between  cyclic  and  monotonic  yield  stress 
fracture  toughness 

crack  tip  stress  intensity  based  on  maximum  cyclic  stress 
effective  value  for 

crack  tip  stress  intensity  based  on  overload  (=  for  CA-loading) 

crack  tip  stress  intensity  based  on  minimum  cyclic  stress 

crack  tip  stress  intensity  based  on  opening  stress  in  cycle 

crack  tip  stress  intensity  based  on  underload  (=  K„|„  for  CA-loading) 

crack  growth  delay  (in  cycles)  after  a  single  overload 

stress  ratio,  such  as  minimum  over  maximum  stress  in  any  cycle,  though 

denoted  as  K^i„  /K„„,o,  in  McAir  model 

size  of  plastic  zone  due  to  overload 

factor  that  accounts  for  material  thickness  in  calculation  of 

factor  that  accounts  for  cyclic  material  behaviour  in  calculation  of 

factor  that  relates  K„in.eff  to  stress  ratio  R 


a  retardation  parameter,  theoretically  representing  material  thickness; 

in  practice  used  to  tunc  the  crack  growth  model 
P  equalto(a+ 1)/2 

AK  stress  intensity  range  (K„„  -  K„i„) 

AK^ff  effccUve  stress  intensity  range  -  K^jn  ^ff) 

maximum  cyclic  stress 

oi  maximum  stress  in  overload  cycle 

Oy  eff  effective  monotonic  yield  stress,  used  in  calculation  of  w^,  to  cover  large 

scale  yielding  of  net  section 
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1.  INTRODUCTION 

The  application  of  crack  growth  predictions  has  been  highly  stimulated  by  the  current 
airworthiness  requirements  imposed  by  civil  and  military  agencies.  These  requirements  were  initiated  to 
ensure  the  fatigue  safety  and  durability  of  a  fleet  of  aircraft  Small  cracks,  either  pre-existing  at  the  time 
of  manufacture  or  created  by  in-service  conditions,  are  assumed  to  grow  during  the  life  of  the  aircraft. 
The  growth  of  the  crack  should  be  predictable  to  provide  guidance  for  inspection  programs,  which 
ensure  that  cracks  will  never  propagate  to  failure  prior  to  detection.  For  this  purpose  many  crack  growth 
models  have  been  devised  during  the  last  three  decades.  None  of  these  models  is  universally  accepted 
however. 

Finney  and  Machin  [1]  have  examined  four  well-known  load-sequence-effect  models,  Wheeler, 
Willenborg,  Modified  Willenborg  and  a  Crack  Closure  model,  for  their  effectiveness  in  predicting 
fatigue  crack  growth.  Predictions  from  the  models  were  compared  against  experimental  crack  growth 
data  obtained  using  several  reconstitutions  of  a  range-pair  counted  flight-by-flight  fighter  load  sequence. 
Although  the  Wheeler  model  gave  reasonable  results,  none  of  the  models  correctly  predicted  the  trends 
in  experimental  crack  growth  among  the  various  reconstituted  sequences. 

This  present  work  describes  the  assessment  of  a  fifth  crack  growth  model,  viz.  the  McDonnell 
Douglas  Contact  Stress  Model  (’McAir  model').  Essentially  this  model  is  based  on  the  phenomenon  of 
crack  closure  at  positive  stress  levels.  Again  the  above  mentioned  variations  on  a  fighter  load  sequence 
have  been  used  to  estimate  the  accuracy  of  this  model.  The  computer  program  that  incorporates  the 
McAir  model  is  fairly  extensive  and  offers  several  stress  intensity  factor  solutions  and  other  options  such 
as:  surface  flaw,  bearing  stress  correction  for  a  crack  at  a  hole  and  cracks  with  a  residual  stress  intensity 
adjustment  [2],  Only  the  basic  option  of  a  through  crack  in  a  centre  cracked  specimen  is  examined  in 
this  report. 


2.  THE  MCAIR  MODEL 

The  McAir  crack  growth  model  is  based  on  crack  closure,  although  it  is  classified  as  a  contact 
stress  model  by  the  authors.  It  has  been  observed  by  Elber  [3]  and  many  others  that  fatigue  cracks 
become  physically  closed  during  unloading  at  positive  stress  levels.  During  subsequent  uploading  the 
fatigue  cracks  are  only  opened  if  the  load  is  sufficiently  high.  This  can  be  explained  with  the  effect  of 
residual  plasticity  that  is  left  in  the  wake  of  the  crack  tip  as  the  crack  grows  through  prior  plastic  zones. 
A  method  has  been  developed  by  Dill  and  Saff  [4,5]  to  determine  the  displacements  and  the  plastic 
deformations  along  the  crack  surfaces.  In  this  way  they  were  able  to  calculate  crack  closure  under 
constant  amplitude  loading  as  a  function  of  the  .stress  ratio.  In  the  McAir  model  this  curve  is 


manipulated  to  obtain  the  minimum  effective  stress  intensity  factor,  or  opening  stress  intensity,  for  each 
block  of  cycles  of  a  variable  amplitude  load  sequence.  It  is  assumed  then  that  the  same  effective  stress 
intensity  range  defined  as  the  maximum  stress  intensity  minus  the  minimum  effective  stress 
intensity  ^  cycle,  under  constant  amplitude  loading  and  under  variable  amplitude  loading  will 

produce  the  same  crack  extension  (similarity  concept). 

2.1  Assumptions 

Like  any  other  crack  closure  model  the  McAir  model  is  based  on  the  assumption  that  crack  growth 
behaviour  is  controlled  by  displacements  within  the  plastic  zone.  These  displacements  are  affected  by 
crack  surface  contact  behind  the  crack  tip  resulting  from  residual  displacements  left  on  the  crack  surface 
as  the  crack  grows  through  prior  plastic  zones.  The  effect  of  crack  surface  contact  during  a  load  cycle  is 
to  restrict  the  change  in  crack  opening  displacements  during  the  cycle,  reducing  the  portion  of  the  load 
cycle  effective  in  propagating  the  crack  [2].  With  this  concept  crack  growth  retardation  after  overloads 
and  acceleration  after  underloads  can  be  explained. 

It  is  assumed  that  for  any  load  cycle,  the  load  level  at  which  closure  occurs  is  a  function  of  the 
highest  peak  load  and  lowest  valley  load  which  have  occurred  prior  to  that  cycle.  This  assumption  only 
holds  if  the  applied  load  sequence  is  more  or  less  random,  however,  with  frequent  occurrence  of 
overloads.  To  provide  for  more-structured  sequences,  the  'highest'  and  'lowest'  load  are  regularly  updated 
as  the  crack  grows  through  the  zone  of  plasticity  caused  by  the  overload.  It  is  also  assumed  that  the 
definition  of  a  load  cycle  is  unimportant,  or,  in  other  words,  that  it  does  not  matter  whether  the  minimum 
cyclic  load  proceeds  or  succeeds  the  maximum  load  in  the  cycle.  Again,  this  only  holds  for  load 
sequences  with  frequent  changes  in  load  amplitude. 

Unlike  most  other  crack  growth  models,  it  is  assumed  that  any  sbess  variation,  no  matter  how 
small,  contributes  to  crack  growth  (provided  that  the  maximum  cyclic  stress  is  greater  than  zero).  In 
other  words,  it  is  assumed  that  some  plastic  deformation  occurs  at  the  crack  lip  even  when  the  crack  is 
only  partially  open.  It  is  left  to  the  user  to  eliminate  insignificant  cycles  from  the  input  load  sequence. 

The  crack  length  that  is  used  for  the  calculation  of  the  stress  intensity  is  assumed  to  be  constant 
during  a  block  of  load  cycles.  This  implies  that  the  number  of  cycles  in  such  a  block  should  be  limited  to 
avoid  inaccurate  answers  in  case  of  highly  loaded  specimens. 

No  distinction  is  made  between  multiple  and  single  overloads  or  underloads,  so  the  effect  of 
multiple  overloads  producing  more  retardation  than  a  single  overload  is  not  accounted  for. 


2.2  Constant  Amplitude  Loading 


The  plastic  deformations  left  in  the  wake  of  tire  crack  tip  cause  the  emek  to  clo.se  before  complete 
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unloading.  Under  constant  amplitude  loading,  the  stress  level  at  which  this  occurs  is  mainly  a  function  of 
the  applied  stress  ratio  R,  though  stress  state  and  material  properties  also  play  a  role.  From  advanced 
analysis  Dill  and  Saff  [4,5]  have  derived  the  locus  of  these  minimum  effective  stress  intensity  factors 
^  example  of  such  a  stress  ratio  correction  curve  is  given  in  figure  2.1.  According  to  Dill  and 
Saff  this  curve  is  lower  for  plane  strain  than  for  plane  stress,  higher  for  cyclic  strain  hardening  materials, 
and  lower  for  strain  softening  materials. 

The  equations  that  describe  the  stress  ratio  correction  curve  for  a  specific  material  have  been 
derived  from  the  listing  of  computer  program  CNTKM8  [6].  This  program  is  the  most  recent  version 
available  to  ARL  of  computer  program  CNTKSPC  as  described  in  reference  [2].  According  to  [6]  the 
minimum  effective  stress  intensity  factor  for  a  stress  ratio  R  that  is  greater  than  or  equal  to  zero 

is  given  by: 

Kmin.eff/K™.x  =  (forR>0)  (2.1a) 

where 

X,  =  0.46733  +  0.29401  R  +  0.23866  R^ 
x„  =  1-(1  Wp)(l.R) 

Xf,  =  1  ■  (0.756  -  0.912fr  +  0.156fr2)(l  -  R) 

P  =  (a+l)/2 

and  is  the  stress  intensity  based  on  peak  cyclic  load.  If  the  suess  ratio  is  less  than  zero  the 
minimum  effective  stress  intensity  can  be  found  with: 


KTnin,eff/K„«  =  efnu-e00S'^  + 0.225  R  (forR<0)  (2.1b) 

where  effkr  is  equal  to  the  ratio  as  given  by  equation  (2.1a)  for  a  stress  ratio  R  that  equals 

zero.  If  equation  (2.1)  yield  a  minimum  effective  .sucss  intensity  that  is  lower  than  zero  or  lower 

than  the  minimum  stress  intensity  in  the  current  cycle,  then  is  set  equal  to  zero  or  K^j„, 

whichever  is  greater. 

Par:.meler  relates  the  minimum  effective  stress  intensity  to  the  stress  ratio  R.  It  is  similar  to 
rclatiors  proposed  by,  or  derived  from,  Elbcr  [3],  Schijve  [7]  and  others. 

Parameter  x„  accounts  for  the  effect  of  material  thickness.  It  depends  on  the  retardation  parameter 
a,  which  is  input  by  the  user  and  ranges  from  -1  (no  retardation)  to  +1  (maximum  retardation). 
Theoretically,  a  value  of  zero  pertains  to  the  situation  of  plane  strain  while  a  value  of  +1  pertains  to  the 
state  of  plane  suess  around  the.  crack  tip.  In  practice,  however,  factor  a  is  used  to  tunc  the  model. 

Parameter  Xf,  brings  into  account  the  behaviour  of  the  material  under  cyclic  loading.  It  is  a  function 
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of  nuilo  fr  between  the  cyclic  end  monotonic  yield  stress.  Cyclic  stnin  hardening  behaviour  will  result  in 
more  retaidatioa  while  cyclic  softening  will  cause  less  retardation.  This  may  be  rather  surprising  since 
retvdation  usoaliy  is  associated  with  the  plastic  tone  size,  which,  in  its  turn,  is  inversely  proportional  to 
the  yield  stress.  It  should  not  be  forgotten  however,  that  the  outer  regions  of  the  plastic  zone  cover  areas 
that  have  hardly  experienced  reversed  plasticity  so  the  plastic  zone  size  is  principally  related  to  the 
monotonic  yield  stress.  Crack  growth,  on  the  other  hand,  is  controlled  by  displacements  within  the 
plastic  zone  rtear  the  crack  tip,  where  the  mate  rial  has  experienced  severe  reversed  plasticity.  So  crack 
growth  will  be  determined  by  the  cyclic  behaviour  of  the  material. 

The  minimum  effective  stress  intensity  K„jn,efr  “  equation  (2.1)  is  depicted  in  figure  2.2 

as  a  function  of  R  and  cl  Figure  2.3  shows  how  K„i„  varies  with  ratio  fr.  Figure  2.4  shows  how  the 
McAir  equations  compare  with  equations  proposed  by  Elbcr  [3]  and  Schijvc  [7]. 

2J  Variable  Amplitude  Loading 

The  stress  ratio  correction  curve  does  not  represent  the  behaviour  at  the  crack  tip  under  variable 
amplitude  loading.  The  minimum  effective  stress  intensity  K^jn^eff  required  to  open  the  crack  is  not  only 
determined  by  the  stress  ratio  of  the  current  cycle,  but  it  also  depends  heavily  on  the  highest  and  lowest 
load  previously  experienced. 

Consider  the  simplified  stress  ratio  correction  curve  of  figure  2.1.  The  same  curve  is  depicted  in 
figure  2.5a  again.  In  this  figure,  however,  this  curve  represents  the  minimum  effective  stress  intensity 
that  would  occur  if  the  material  were  subject  to  a  constant  amplitude  loading  with  maximum  stress 
intensity  and  minimum  stress  intensity  where  and  are  defined  by  a 

characteristic  ovcrload/underload  combination  prior  to  the  current  cycle  of  the  actual  sequence.  The 
definition  of  these  loads  will  be  dealt  with  later  on  in  this  chapter.  To  find  the  minimum  effective  stress 
intensity  of  the  current  load  cycle,  point  A  on  the  stress  ratio  correction  curve  that  is  associated 

with  the  ovcrload/underload  combination  is  used  as  a  starting  point.  found  by  drawing  a  line 

with  tangent  0.225  to  point  B  that  corresponds  with  the  ratio  /K^„  p,  on  the  abscis,sa,  where 
is  the  minimum  applied  stress  intensity  in  the  current  cycle.  It  may  well  happen  that  the  resulting 
minimum  effective  stress  intensity  K^u,  is  lower  than  zero  or  lower  than  the  minimum  applied  suess 
intensity  K^,„,  a.s  is  shown  in  figure  2.5c.  In  that  case  eff  is  assumed  to  be  equal  to  zero  or  to 
whichever  is  greater.  It  should  be  noted  that  the  ordinate  represents  the  ratio  of  and  q,,  and 

thus  It  is  clear  that  the  minimum  effective  stress  intensity  K^i„efr  in  the  current  cycle  depends  heavily 
on  the  overload  K„„  that  occurred  before.  A  high  overload  will  give  a  high  value  for  If  the 

tangent  of  line  AB  were  zero  instead  of  0.225  then  would  .solely  be  determined  by  the 

ovcrload/underload  combination. 

The  effect  of  the  undcrlrjad  is  Illustrated  with  figure  2.5b.  If  the  underload  goc,s  down  from 
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to  K„i„,or>  keeping  K„jn  and  K„„,oi  “me  level,  then  the  minimum  effective  stress  intensity 
^in,eff  80CS  down  from  B  to  B'.  This  reflects  the  effect  of  underloads  to  reduce  crack  closure.  It  is 
interesting  to  note  that  the  tangent  of  the  stress  ratio  curve  for  R  <  0  is  approximately  equal  to  0.225,  as 
can  be  seen  in  figure  2.2.  This  means  that  the  actual  magnitude  of  negative  underloads  does  not  really 
matter  for  crack  growth  acceleration. 

A  peculiar  feature  of  the  McAir  model  is  the  introduction  of  a  maximum  effective  stress  intensity 
analogous  to  the  minimum  effective  stress  intensity  Usually  K^^eff  ** 

maximum  stress  intensity  that  occurs  in  a  cycle,  as  one  would  expect  If  the  overload  is  much  higher 
than  this  maximum  cyclic  stress,  however,  calculated  such  that  the  effective  stress  intensity 

range  AK^ff  in  the  cycle  becomes  equal  to  0.225  (K,nj„  -  K^bi).  where  and  K^i„  are  the  actual 
cyclic  stresses.  In  this  way  sustained  crack  arrest  due  to  high  overloads  never  occurs. 

2.3.1  Definition  of  Ovaload  and  Underload 

As  already  outlined  the  closure  stress  for  a  particular  load  cycle  is  mainly  determined  by  a 
characteristic  overload  and  underload  which  have  occurred  prior  to  that  cycle.  These  loads  are  regularly 
updated  as  the  crack  grows  through  the  plastic  zone  caused  by  the  overload.  The  overload  is  adjusted  to 
the  maximum  load  in  the  current  cycle  if  (1)  the  crack  has  grown  through  the  plastic  zone  caused  by  the 
previous  overload,  or  (2)  the  maximum  stress  intensity  in  a  cycle  exceeds  the  value  of  a  variable  named 
check.  This  variable  is  a  function  of  the  previous  overload,  the  size  of  the  associated  plastic  zone  w^,  and 
crack  length  a^i,  and  the  current  crack  length  a.  This  is  shown  in  figure  2.6,  which  was  derived  from 
reference  [6].  Also  indicated  in  this  figure,  for  comparison,  is  the  stress  intensity  factor  necessary  to 
create  a  plastic  zone  that  extends  beyond  the  plastic  zone  as  created  by  the  previous  overload.  This  curve 
was  derived  by  assuming  that  the  size  of  the  plastic  zone  due  to  some  arbitrary  load  is  proportional  to 
K^,  where  K  is  the  stress  intensity  associated  with  this  load  (see  equation  (2.2)  in  the  next  section). 

In  both  above-mentioned  cases  the  underload  is  updated  as  well.  Apparently  the  physical  reasoning 
behind  this  is  that  an  overload  erases  the  material's  memory.  The  underload  is  also  adjusted  to  the 
minimum  load  in  the  current  cycle  when  the  minimum  stress  intensity  in  a  cycle  is  lower  than  the  stress 
intensity  that  is  associated  with  the  previous  underload.  In  short; 


is  adjusted  to  when 

1. 

a  >  aol  +  '^oi 

2. 

>  check 

is  adjusted  to  when 

1. 

a  >  ^oi  +  w„, 

2. 

>  check 

3. 

The  above  set  of  criteria  is  summarized  in  figure  2.7,  which  shows  a  very  simple  flow  diagram  of 
the  McAir  model.  The  notation  is  the  same  as  used  throughout  the  computer  program,  and  is  explained 
in  appendix  A.  Again,  only  the  case  of  a  through  crack  in  a  CCT-specimen  is  considered. 

2.3.2  Thickness  Effect 

In  theory  the  retardation  parameter  a  accounts  for  the  stress  state  iL-ound  the  crack  tip.  A  value  of 
+1  would  pertain  to  thin  sheet  material  and  a  value  of  zero  to  thick  specimens.  Factor  a  asserts  its 
influence  in  two  ways.  First  of  all  a  determines  the  locus  of  the  stress  ratio  correction  curve,  as  is  clear 
from  equation  (2.1)  and  from  figure  2.2.  This  means  that  a  affects  crack  growth  under  both  constant  and 
variable  amplitude  loading.  Secondly  factor  a.  determines  the  size  of  the  plastic  zone  ahead  of  the  crack 
tip.  This  directly  affects  crack  growth  after  an  overload,  since  the  effect  of  an  overload  does  not  vanish 
until  the  crack  has  grown  through  the  associated  plastic  zone. 

To  determine  the  importance  of  the  stress  ratio  correction  curve,  we  will  first  consider  the  case  of 
no-load-interaction  during  spectrum  crack  growth.  In  that  case  the  stress  ratio  correction  curve  is  solely 
used  to  account  for  the  stress  ratio  effect  on  the  opening  stress.  The  effect  of  a  on  crack  growth  need  not 
necessarily  be  very  large  then,  although  the  range  of  curves  in  figure  2.2  suggests  otherwise.  This  has  to 
do  with  the  way  the  McAir  model  treats  the  material  crack  growth  curve.  Starting  point  is  the  da/dN-AK 
curve  for  a  stress  ratio  R  of  zero,  which  must  be  input  by  the  user  of  the  model.  The  tabulated  AK- 
values  are  then  transformed  to  AK^ff- values  with  the  help  of  equation  (2.1),  where  the  effective  stress 
intensity  range  AK^jf  is  equal  to  minus  It  is  obvious  from  figure  2.2  that  lower  AK^fp 

values  will  be  obtained  for  thin  material  (a  =  1)  than  for  thick  material  (a  =  0).  During  the  subsequent 
crack  growth  calculation  AK^jj  is  calculated  for  each  block  of  cycles.  Again  it  is  obvious  that  lower 
values  for  AK^fj  will  be  obtained  for  thin  material  than  for  thick  material.  The  associated  da/dN  -  value 
then  is  found  by  interpolation  in  the  previously  obtained  da/dN  -  AK^ff  table.  Because  factor  a  affects 
both  the  input  crack  growth  data  and  the  stress  intensity  range  due  to  an  applied  load  cycle,  the  resulting 
effect  on  stress  ratio  correction  is  quite  small.  This  is  illustrated  with  the  following  example.  Consider 
the  fictitious  though  realistic  da/dN  -  AK  data  set  as  listed  in  the  table  below: 


da/dN 

[m/cycle] 

[MPaVm] 

a  =  0 

a  =  1 

110-’ 

3-10-® 


3.0 

30.0 


2.0 

20.0 


1.6 

16.0 
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The  flm  two  columns  represent  the  input  date  for  e  stress  ratio  of  zero  while  the  last  two  columns 
represent  the  anocieted  values  for  two  extremes  of  factor  a  as  calculated  with  equation  (2.1a), 
with  retk)  fr  equal  to  one.  Now  consider  an  arbitrary  load  cycle  with  a  stress  ratio  R  of  0.5  and  a  stress 
intensily  range  AK  of  IS.O  MPaVm.  Application  of  equation  (2.1a)  to  this  value  yields  an  effective  stress 
intensity  range  of  12.74  MPaVm  for  a  equal  to  zero  and  9.78  MPaVm  for  a  equal  to  one 
respectively.  Interpolation  in  the  respective  AK^«  columns  of  the  above  uble,  using  the  logarithm  of 
the  listed  values  as  is  usual,  finally  yields  a  growth  rate  da/dN  of  6.3*10''’  m/cycle  for  a  equal  to  zero  and 
m/cycle  for  a  equal  to  one  respectively.  The  difference  between  the  calculated  growth  rates 
amounts  to  13%,  which  is  relatively  small  considering  the  fairly  large  spread  between  the  stress  ratio 
curves  for  the  two  extremes  of  factor  a  (see  figure  2.2).  It  can  thus  be  concluded  that  it  docs  not  matter 
much  what  value  for  a  is  used  when  correlating  constant  amplitude  crack  growth  under  arbitrary  stress 
ratio  to  crack  growth  data  for  a  stress  ratio  of  zero,  as  long  as  factor  a  is  applied  consistently.  If  this  is 
not  the  case  large  errors  should  be  expected. 

An  interesting  consequence  of  the  large  spread  between  the  various  stress  ratio  correction  curves  is 
that  the  McAir  model  predicts  different  constant  amplitude  crack  growth  behaviour  for  thin  specimens 
and  thick  specimens.  Consider  again  the  material  crack  growth  curve  for  a  stress  ratio  of  zero  that  is 
given  by  the  table  above.  Suppose  this  curve  has  been  derived  from  tests  on  a  thin  specimen,  so  that  the 
effective  stress  intensity  range  is  given  by  the  last  column  of  this  table.  Now  let  us  apply  this  crack 
growth  curve  to  a  thick  specimen  that  is  subjected  to  constant  amplitude  loading  with  a  stress  ratio  equal 
to  zero  and  a  stress  intensity  range  AK  of  3.0  MPa'^m.  The  effective  stress  intensity  range  AK^jj  then  can 
be  found  with  equation  (2.1a),  using  a  value  for  a  of  zero,  and  is  equal  to  2.0  MPa-^m.  Interpolation, 
using  the  logarithm  of  the  listed  values  of  the  last  and  first  column  respectively,  yields  a  crack  growth 
rate  da/dN  of  2.2-10'’,  while  a  value  of  I.O'IO'*  would  have  been  obtained  for  a  thin  specimen.  It  is  not 
clear  from  publications  whether  such  a  large  difference  is  realistic  or  not,  but  it  is  clear  that  care  should 
be  taken  when  choosing  a  crack  growth  curve  for  a  specific  problem. 

Regarding  variable  amplitude  loading,  both  the  locus  of  the  sucss  ratio  correction  curve  and  the 
size  of  the  plastic  zone  arc  of  importance.  The  equation  that  is  used  for  the  plastic  zone  size  w^,  due  to 
an  overload  K„„  „,  is: 


where 


(2.2) 


p  =  (a+  l)/2 

effective  monotonic  yield  stress  by  which  large  scale  net  section  yielding  is  accounted 
for.  At  lower  stres*;  levels  is  approximately  equal  to  the  monotonic  yield  .stress,  while  its  saluc 
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drops  drastically  when  the  applied  stress  approaches  the  monotonic  yield  stress.  This  is  shown  in  figure 
2.8.  Equation  (2.2)  yields  a  plastic  zone  size  (hat  is  twice  as  large  under  plane  strain  as  under  plane 
stress.  Usually  this  ratio  is  assumed  to  be  equal  to  three  [8].  It  was  already  discussed  in  section  2.3.1  that 
in  the  case  of  a  load  sequence  with  infirequent  changes  in  load  amplitude  the  effect  of  an  overload  does 
not  vanish  until  the  crack  has  grown  through  the  associated  plastic  zone.  Thus  the  thickness  of  a 
specimen  will  play  an  important  role  in  the  process  of  crack  growth  retardation  after  an  overload.  This 
effect  is  promoted  by  the  dependence  of  the  minimum  effective  stress  intensity  in  an  overload 

cycle  (which  is  given  by  the  stress  ratio  correction  curve)  on  parameter  a.  The  higher  the  value  of  a,  or 
the  smaller  the  specimen  thickness,  the  higher  the  minimum  effective  stress  intensity  in  the  overload 
cycle,  and  thus  the  smaller  the  effective  stress  intensity  range  in  subsequent  cycles. 

2.3.3  Examples 

Two  examples  of  the  definition  of  overloads  and  underloads  are  given  in  figure  2.9.  The 
consecutive  overloads  and  underloads  are  indicated  with  arrows.  The  approximate  location  of  the 
minimum  effective  stress  intensity  is  indicated  with  a  small  cross.  The  initial  overload/underload 
combination  is  determined  by  the  first  cycle,  which  is  sclf-cvidenu  According  to  the  above  criteria  the 
second  cycle  in  figure  2.9a  yields  a  new  overload,  since  the  maximum  stress  intensity  associated  with 
this  cycle  will  exceed  the  value  of  variable  check  (see  figure  2.6).  As  a  consequence  a  new  underload  is 
defined  as  well,  since  this  is  always  coupled  to  the  definition  of  a  new  overload.  At  the  third  cycle 
nothing  happens.  Even  the  minimum  effective  stress  intensity  remains  the  same,  since  equals 
^in.oi-  fourth  cycle  yields  a  new  underload  and  the  minimum  effective  stress  intensity  goes  down 
as  a  consequence.  The  fifth  cycle  is  dominated  by  the  combination  of  the  cycle  2  overload  and  the  cycle 
4  underload.  The  minimum  effective  stress  intensity  is  higher  than  that  of  the  fourth  cycle,  but  it  will  not 
become  as  high  as  the  minimum  effective  stress  intensity  of  the  third  cycle. 

The  coupling  of  the  definition  of  underloads  to  the  definition  of  overloads  can  be  rather  loose  in 
certain  situations.  This  is  illustrated  in  figure  2.9b.  The  effect  of  the  severe  underload  at  the  second  cycle 
is  immediately  cancelled  out  by  the  overload  at  the  third  cycle,  which  is  not  very  realistic. 

It  should  be  noted  that  the  McAir  model  does  not  explicitly  bring  into  account  the  effect  of 
multiple  overloads  producing  more  retardation  than  a  single  overload.  This  effect  can  only  be  arbitrarily 
accounted  for  by  tuning  the  model  with  retardation  parameter  a.  It  should  also  be  noted  that  the  effect  of 
delayed  retardation  after  a  high  overload  is  not  included  in  the  model. 
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3.  ASSESSMENT  OF  THE  MCAER  MODEL 


3.1  Experiments 

To  assess  the  McAir  model  use  has  been  made  of  experimental  data  described  in  references  [1]  and 
[9].  These  data  were  obtained  from  tests  using  several  reconstitutions  of  a  flight-by-flight  fighter  load 
sequence.  This  sequence  was  that  used  in  fatigue  testing  a  Mirage  IIIO  wing.  It  is  typical  of  that  for  a 
fighter  aircraft,  covering  500  flights  and  consisting  of  108,920  turning  points.  The  maximum  load 
corresponds  to  7.5  g.  Cycles  in  this  sequence  were  counted  using  the  range-pair  method,  giving  a  table 
of  numbers  corresponding  to  various  peak-trough  combinations.  For  the  present  work  the  original  flight- 
by-flight  sequence  plus  three  artificial  reconstitutions  of  this  sequence  were  used:  Two  more-structured 
sequences  were  supposed  to  promote  and  retard  crack  growth.  The  acceleration  sequence  aimed  at 
shortening  of  crack  growth  life,  while  the  retardation  sequence  was  developed  to  produce  the  longest 
possible  crack  growth  life.  These  two  sequences  are  shown  in  figure  3.1.  The  fourth  reconstitution,  the 
random  sequence,  was  generated  in  a  more  or  less  random  way.  All  sequences  used  returned  the  same 
range-pair  table  when  counted.  Information  about  the  way  these  load  sequences  were  generated  can  be 
found  in  reference  [1]. 

Centfe-cracked  specimens  of  A7-U4SG-T651  aluminium  alloy  (equivalent  to  US  alloy  2214)  were 
tested  under  each  of  the  four  sequences  described  above.  The  specimens  measured  12  mm  x  80  mm  x 
300  mm,  with  a  central  spark-machined  slit  0.015  mm  wide  by  12  mm  in  length.  The  specimens  were 
pre-cyclcd  to  a  crack  length  of  13.2  mm.  The  average  lest  frequency  was  5.5  Hz,  and  two  gross-area 
stress  levels  of  11.09  and  17.(K)  MPa/g  were  employed  to  yield  a  range  of  test  lives.  More  details  about 
the  specimens  and  the  test  procedures  can  be  found  in  references  [9]  and  [10]. 

The  lest  results  are  presented  in  figure  3.2,  which  has  been  obtained  from  reference  [1],  A 
remarkable  result  is  that  the  acceleration  sequence  yielded  the  longest  crack  growth  life  under  both  stress 
levels.  Furthermore,  the  differences  in  crack  growth  lives  between  those  under  the  acceleration  and  the 
retardation  sequences  and  also  between  those  under  the  random  and  f-b-f  sequences  were  quite  small. 
Crack  growth  under  the  more-structured  sequences  (acceleration  and  retardation),  however,  is 
significantly  slower  than  under  the  less-structured  more-fluctuating  sequences  (random  and  f-b-f).  The 
maximum  difference  in  average  crack  growth  lives  was  between  those  under  the  acceleration  and 
random  sequences  and  was  an  average  of  1.59:1  at  the  stress  level  of  11.09  MPa/g,  and  1.51:1  at  17.0 
MPa/g. 


3.2  Constant  Amplitude  Data 


The  McAir  model  assumes  that  the  constant  amplitude  da/dN- AK  curve  that  is  input  by  the  u.scr  is 
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for  a  stress  ratio  of  zero.  Only  experimental  data  for  stress  ratio's  -0.7,  0.2  and  0.7  were  available  for  the 
A7-U4SG-T651  aluminium  alloy  however  [11].  To  overcome  this  problem  an  equation  due  to  Schijve 
[7]  was  used  to  define  an  effective  AK  in  such  a  way  that  it  correlates  all  the  experimental  data 
available.  This  equation  had  already  been  applied  successfully  in  the  past  at  ARL.  It  can  be  written  as; 

AKrff  =  (0.55  +  0.35R  +  0.1  R^)  AK  (3.1) 

This  equation  is  represented  in  figure  2.4  by  the  'Schijve  (1)'- curve.  Essentially  this  equation  is  similar 
to  equation  (2.1).  Equation  (3.1)  showed  better  correlation  between  the  constant  amplitude  data, 
however.  With  equation  (3.1)  the  three  sets  of  experimental  data  were  transformed  to  one  set  of  R  =  0 
data.  Using  logarithmic  scales  a  straight  line  was  fitted  through  these  data.  This  resulted  in  the  following 
material  curve,  which  is  valid  for  crack  growth  rate  in  m/cycle  and  stress  intensity  in  PaVm: 

log(da/dN)  =  3.7162  log(AK)- 33.158  (forR  =  0)  (3.2) 

Two  curves  were  added  to  this  straight  line  to  account  for  the  growth  rates  outside  the  AK  data  limits. 
These  curves  were  fitted  in  a  rather  arbitrary  way,  taking  care  however  that  the  tangent  of  these  curves  at 
the  AK  data  limits  was  equal  to  the  tangent  of  the  straight  line  through  the  data.  The  threshold  value  of 
the  effective  AK  was  taken  to  be  0.80  MPaVm  with  da/dN  equal  to  lO*”  m/cycle,  and  the  fracture 
toughness  was  taken  as  48.6  MPaVm  with  da/dN  equal  to  1.0  m/cycle.  These  values  were  obtained  from 
reference  [10].  The  resulting  crack  growth  curve  is  shown  in  figure  3.3.  The  tabulated  data  are  given  in 
table  1.  Table  2  gives  the  static  properties  of  the  material  used.  These  data  were  obtained  from 
references  [10]  and  [12]. 

3.3  Predictions 

Crack  growth  was  predicted  using  the  McAir  model  and  the  four  sequences  described  in  section 
3.1,  for  the  two  stress  scaling  levels  of  11.09  and  17.0  MPa/g.  The  da/dN-AK  data  of  table  1  and  the 
static  properties  of  table  2  were  used.  It  was  not  clear  what  value  for  the  input  parameter  a  should  be 
used.  Given  the  relatively  large  specimen  thickness  of  12  mm,  a  value  of  zero  seemed  appropriate  since 
this  value  theoretically  pertains  to  the  situation  of  plane  strain  around  the  crack  lip.  Using  various  values 
for  a  showed  that  much  lower  values  yielded  belter  predictions,  especially  for  the  stress  scaling  level  of 
1 1 .09  MPa/g.  It  was  fell  that  nothing  is  proved  by  getting  excellent  predictions  due  to  an  input  parameter 
that  wa.s  manipulated  to  suit  experimental  rc.sulLs.  Therefore  it  was  decided,  rather  arbitrarily,  to  use  an 
a  value  of  -0.2  for  both  stress  levels,  as  a  compromi.se  between  experiment  and  theory.  The  computer 
program  CNT8SHR,  a  modified  version  of  the  originally  available  computer  program  CNTKM8  [6],  was 
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used  for  the  calculations.  This  program  is  described  in  appendix  A. 

Predicted  crack  growths  arc  shown  in  figure  3.4  for  the  four  sequences  and  the  two  stress  levels 
examined.  In  figure  3.5  the  predictions  from  the  McAir  model  are  compared  against  the  experimental 
crack  growth  data.  This  figure  also  shows  the  predictions  for  the  unretarded  or  no-load-interaction-casc. 
These  latter  predictions  were  made  using  a  simple  computer  program  that  employed  equation  (3.1)  to 
account  for  the  stress  ratio  effect  in  each  cycle.  It  would  also  have  been  possible  to  use  computer 
program  CNT8SHR,  or  CNTKM8,  to  make  predictions  for  the  unretarded  case  by  inputting  a  value  of  -1 
for  a.  This  would  not  have  yielded  the  same  results  however,  since  no  correction  for  the  stress  ratio  is 
made  if  a  equals  -1  and  all  cycles  are  treated  as  if  the  stress  ratio  were  zero. 

Table  3  summarizes  the  predicted  crack  growth  lives  for  various  values  of  retardation  parameter  a. 
The  predictions  for  the  unretarded  case  are  also  listed.  Note  the  difference  between  these  latter 
predictions  and  the  predictions  with  a  equal  to  -1.  Table  4  essentially  contains  the  same  information  as 
table  3,  except  that  the  predicted  crack  growth  lives  are  presented  as  a  percentage  of  the  experimental 
crack  growth  life.  In  table  5  the  predicted  crack  growth  lives  are  compared  against  predictions  made 
with  three  other  crack  growth  models:  the  Wheeler  model,  the  Willenborg  model  and  the  Modified 
Willenborg  model.  The  numbers  between  brackets  denote  the  predicted  crack  growth  lives  as  a 
percentage  of  the  experimental  life.  These  data  were  obtained  from  reference  [1].  The  data  arc  listed  in 
the  order  of  increasing  predicted  crack  growth  lives. 

Several  observations  follow  from  the  predicted  data.  Table  4  shows  that  predictions  for  the 
unretarded  case  were  always  less  than  experimental  lives,  the  difference  for  the  more-random  sequences 
at  the  lower  stress  level  being  quite  small,  however.  For  these  sequences  at  this  sucss  level  die 
predictions  with  the  McAir  model,  with  a  equal  to  -0.2,  yielded  poor  results.  The  crack  growth  lives  for 
these  cases  were  grossly  over-estimated  by  a  factor  of  two.  An  interesting  detail  is  that  even  with  a  value 
for  a  as  low  as  -1  the  crack  growth  life  is  still  over-estimated.  The  crack  growth  lives  for  the  more- 
suucturcd  sequences  at  the  lower  stress  level  were  predicted  very  well.  The  maximum  deviation  from 
the  experimental  results  was  1 1  %.  The  predictions  over  the  whole  range  of  sequences  at  the  higher  stress 
level  were  fairly  good.  The  maximum  difference  between  predicted  and  experimental  lives,  using  a 
proportion  greater  than  one,  is  1.39:1  for  the  acceleration  .sequence.  It  should  be  realized  that  these 
results  arc  highly  affected  by  the  value  for  input  parameter  a.  Apparently  the  value  of  -0.2  is  very 
satisfactory  for  the  stress  scaling  level  of  17.0  MPa/g.  For  the  stress  level  of  1 1.09  MPa/g  a  lower  value 
for  a  would  be  more  appropriate  to  get  a  better  average  result  over  the  whole  range  of  load  sequences. 
For  reasons  mentioned  before  no  attempt  has  been  made  to  optimize  a  for  this  stress  level  however.  In 
that  respect  the  ranking  of  the  predictions  with  the  McAir  model  in  table  5  should  not  be  taken  too 
absolutely. 

Obtaining  reasonable  crack  growth  predictions  may  be  important  for  normal  engineering  purposes, 
but  whether  or  not  the  right  trend  is  predicted  is  of  much  more  importance  when  as.scssing  a  crack 
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growth  model.  In  that  respect  it  is  interesting  to  compare  the  various  predictions  among  themselves. 
Figure  3.4  shows  that  a  longer  crack  growth  life  was  predicted  for  the  acceleration  sequence  than  for  the 
retardation  sequence,  though  the  difference  is  quite  small.  This  is  in  agreement  with  the  experimental 
data,  as  can  be  appreciated  from  figure  3.2.  Comparison  between  the  results  for  the  more-random 
sequences  also  shows  a  good  correlation  between  experimental  and  predicted  trend.  However, 
experimentally  the  acceleration  sequence  gave  the  longest  crack  growth  life  and  the  random  sequences 
gave  the  shortest  Yet  the  McAir  model  predicted  the  reverse  at  both  stress  levels,  as  can  be  seen  in 
figure  3.4.  According  to  reference  [1]  the  Wheeler  model,  the  Willenborg  model  and  the  Modified 
Willenborg  model  showed  the  same  behaviour. 

3.4  Thickness  Effect 

It  has  been  reported  by  various  authors  that  the  specimen  thickness  seriously  affects  the  fatigue 
crack  growth  behaviour  under  variable  amplitude  loading,  crack  growth  lives  for  thicker  material  being 
systematically  shorter.  Schijve  [13]  has  given  an  overview  of  published  results  of  flight  simulation  tests. 
He  found  differences  of  up  to  a  factor  6.2  between  the  crack  growth  lives  of  thick  and  thin  specimens, 
though  the  variation  of  this  factor  was  large.  Revill  et  al.  [14]  tested  the  effect  of  a  single  overload  on 
subsequent  crack  growth  behaviour  in  2024-T3  aluminium.  They  found  that  the  resulting  delay  in  crack 
growth  Nj  was  greater  by  a  factor  of  8.5  for  1.6  mm  thick  specimens  than  for  6.4  mm  thick  specimens. 
The  magnitude  of  Nj  depends  on  various  factors  like  the  magnitude  of  overload  and  subsequent  cyclic 
loading,  as  can  be  appreciated  from  figure  3.6.  This  figure  is  due  to  De  Koning  and  was  obtained  from 
reference  [15].  It  is  clear  that  the  delay  of  crack  growth  becomes  significantly  larger  with  increasing 
overload. 

To  check  whether  the  McAir  model  correctly  predicts  the  trend  in  experimental  crack  growth  of 
thin  specimens  giving  a  longer  delay  after  a  single  overload  than  thick  specimens,  a  systematic  analysis 
has  been  performed  involving  variation  of  retardation  parameter  a,  representing  specimen  thickness,  and 
overload  attempt  has  been  made  to  reproduce  the  experimental  data  as  published  in 

references  [14]  and  [15].  Only  the  simple  case  of  a  through  crack  in  an  infinite  sheet  under  R  =  0 
loading  has  been  considered.  The  material  data  as  described  in  section  3.2  have  been  used.  The  crack 
dimension  a^i  at  the  time  of  the  overload  was  12.73  mm,  and  the  maximum  stress  in  the  cycles 
following  the  overload  was  50.0  MPa.  Factor  a  ranged  from  -0.4  to  +1.0,  with  an  increment  of  0.2,  and 
the  overload  ratio  ranged  from  1 .05  to  2,95,  where  is  the  maximum  stress  in  the 

overload  cycle. 

The  results  of  the  calculations  arc  presented  in  figure  3.7.  Eight  curves  arc  shown  in  this  figure, 
each  representing  the  crack  growth  delay  period  Nj  as  a  function  of  the  overload  ratio  for  a  specific 
value  of  a.  Several  observations  follow  from  these  curves.  The  most  striking  feature  of  the  curves  is 
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their  shape.  The  first  part  of  each  curve  resembles  a  parabola,  just  like  the  curve  of  figure  3.6.  Each 
curve,  however,  sharply  bends  down  at  a  certain  overload  ratio  to  continue  in  a  much  suaighter  way. 
This  even  applies  to  the  curve  that  is  associated  with  a = -0.4,  though  for  this  curve  the  kink  falls  outside 
the  figure.  The  explanation  is  fairly  simple  and  lies  in  the  definition  of  the  maximum  effective  stress 
intensity  This  stress  intensity  is  normally  equal  to  K,nM*  increasing  overload  the  effective 

stress  intensity  range  (=Kiiux.eff'*^»in.efr)  ^  ^  subsequent  cycles  will  decrease,  and  at  some 
point  it  will  happen  that  in  the  subsequent  cycles  becomes  lower  than  0.225AK.  In  that  case 
is  calculated  such  that  AK^  becomes  equal  to  0.225 AK,  irrespective  of  the  magnitude  of  the 
overload.  In  this  way  sustained  crack  arrest  is  avoided  and  the  crack  will  always  grow  through  the  plastic 
zone  due  to  the  overload.  Since  for  a  certain  overload  ratio  increases  with  increasing  a,  as  can 

be  deduced  from  figure  2.2  and  equation  (2.1),  the  kink  in  the  curve  will  occur  at  lower  overload  levels 
for  higher  values  of  factor  a. 

It  is  not  clear  whether  the  second  part  of  each  curve  represents  the  actual  behaviour  of  a  fatigue 
crack  after  a  high  overload.  Evidence  exists  [16]  that  sustained  crack  arrest  never  occurs  and  that  a 
fatigue  crack  always  will  grow  further,  however  slow.  This  may  explain  why  the  authors  of  the  McAir 
model  introduced  a  minimum  effective  stress  intensity  range,  though  the  value  of  0.225AK  seems 
somewhat  artificial  and  may  in  practice  depend  on  the  specimen  thickness  as  well. 

If  we  disregard  the  effect  of  the  minimum  effective  stress  intensity  range  and  only  look  at  the  first 
part  of  the  curves,  it  is  immediately  clear  that  a  thinner  specimen  (higher  a)  will  show  a  much  longer 
delay  period  after  an  overload  than  a  thick  specimen.  Comparing  the  curves  for  a  =  0  and  a  =  1  it  is 
found  that  the  difference  between  the  predictions  ranges  from  3:1  for  an  overload  ratio  of  1.3  to  9:1  for 
an  overload  ratio  of  1.7.  This  seems  in  agreement  with  the  experimental  results  quoted  before. 


3.5  Discussion 

The  crack  growth  predictions  for  the  structured  sequences  were  quite  satisfactory.  A  reasonable 
value  for  the  retardation  parameter  a  (that  is,  a  value  close  to  zero)  yielded  reasonable  predictions.  This 
is  not  true  for  the  random  sequences  however.  For  both  stfcss  levels  the  McAir  model  over-estimated  the 
crack  growth  lives.  Even  a  value  of  -1  for  a  could  not  remedy  this  for  the  lower  stress  level.  In  fact,  the 
experimental  crack  growth  lives  under  the  random  load  sequences  at  the  stfess  level  of  11.09  MPa/g 
were  quite  close  to  the  no-load-interaction  predictions,  as  can  be  appreciated  from  figure  3.5a.  This 
effect  was  already  observed  by  Schijve  [17].  He  found  that  load  sequence  effects  only  occur  when 
changes  in  load  amplitude  arc  infrequent.  In  case  of  frequent  load  changes,  the  effect  of  overloads 
apparently  is  cancelled  out  immediately  by  frequently  occurring  underloads.  This  raises  the  question 
whether  the  effect  of  underloads  is  sufficiently  being  accounted  for  in  the  McAir  model. 

In  the  previous  section  the  physical  significance  of  the  retardation  parameter  a  is  emphasised.  The 


14 


lower  the  value  for  a,  or  the  larger  the  specimen  thickness,  the  shorter  the  crack  growth  delay  after  a 
single  overload.  This  is  entirely  in  accordance  with  experimental  results.  On  the  other  hand,  the  results 
of  the  crack  growth  predictions  suggest  that  the  requited  value  for  a  also  depends  on  the  stress  scaling 
level.  Higher  stress  levels  tend  to  require  a  higher  value  f<»  this  parameter.  This  illustrates  the  fact  that  it 
is  necessary  to  use  this  factor  as  a  tuning  parameter  in  practice. 

The  assumption  of  the  crack  being  stationary  during  a  block  of  cycles  has  little  effect  on  predicted 
crack  growth  life.  As  an  example  the  crack  growth  life  under  the  retardation  sequence  was  calculated 
cycle-by-cycle  for  the  no-load-interaction  case.  This  sequence  was  the  one  with  the  fewest,  and  thus 
largest,  blocks.  The  result  is  presented  as  a  footnote  in  table  3.  Compared  with  the  'normal'  prediction 
the  difference  only  amounted  3  %.  Apparently  the  number  of  cycles  in  the  blocks  was  still  small  enough 
to  give  sufficiently  accurate  results. 


4.  ENHANCED  CRACK  CLOSURE  MODEL 

Although  the  McAir  model  predicted  the  crack  growth  lives  under  the  structured  sequences  quite 
well,  the  results  for  the  random  sequences  were  fairly  disappointing.  Even  more  bothering  was  the  fact 
that  the  experimental  trend  of  the  structured  sequences  giving  the  longest  crack  growth  lives  was  not 
predicted  correctly.  It  was  felt  that  this  was  caused  by  the  twofold  combination  of  poor  modelling  of  the 
underloads  and  inaccurate  calculation  of  the  effective  stress  intensity  range  during  an  overload  cycle. 
This  will  be  explained  in  the  next  paragraphs.  To  check  this  hypothesis  a  new  model,  the  'ARL  model', 
was  devised  which  finally  resulted  in  computer  program  ARLECC  (ARL  Enhanced  Crack  Closure 
model).  This  present  chapter  describes  the  actual  model,  while  the  computer  program  is  described  and 
listed  in  appendix  B. 

4.1  Actual  Model 

In  conception  the  ARL  model  is  very  similar  to  the  McAir  model.  The  definition  of  overloads  and 
underloads  is  nearly  tlie  same,  and  the  handling  of  the  stress  ratio  correction  curve  to  obtain  the 
minimum  effective  stress  intensity  in  a  cycle  is  practically  equivalent.  The  retardation  parameter  a  has 
also  been  maintained,  and  is  to  be  used  in  the  same  way.  The  main  difference  lies  in  the  handling  of  the 
load  sequence,  the  ARL  model  incorporating  a  cycle-by-cycle  calculation,  and  the  calculation  of  the 
minimum  effective  stress  intensity  in  a  cycle.  Furthermore  a  number  of  details  has  changed,  such  as  the 
calculation  of  the  plastic  zone  size  under  plane  strain. 
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4.1.1  Constant  Amplitude  Loading 

In  principle,  equation  (2.1)  can  be  used  to  correlate  constant  amplitude  crack  growth  data  for 
various  values  of  stress  ratio  R.  Applying  these  equations  to  the  available  da/dN  -  AK  curves  for  the 
A7-U4SG-T651  aluminium  alloy,  however,  yielded  fairly  poor  results.  This  was  already  discussed  in 
section  3.2.  For  this  reason  it  was  decided  to  employ  another  equation  for  the  stress  ratio  correction 
curve.  Two  promising  alternatives  were  available,  both  due  to  Schijve  [7].  They  are  depicted  in 
figure  2.4.  These  equations  had  already  been  shown  to  correlate  experimental  data  very  well.  According 
to  reference  [18],  the  first  one,  denoted  with  'Schijve  (1)'  and  given  by  equation  (3.1),  gave  the  best 
results  when  applied  to  a  particular  set  of  constant  amplitude  data.  The  difference  between  the  two 
equations  was  minimal,  however.  Since  the  second  equation  shows  a  more  pronounced  behaviour  for 
negative  R-values,  as  can  be  appreciated  from  figure  2.4,  it  was  decided  to  employ  this  equation  for  the 
ARL  model.  It  can  be  written  as; 

AKjff  =  (0.55  +  0.33R  +  0.12R2)  (4.1) 

Strictly  speaking  this  equation  is  only  valid  for  2024 -T3  aluminium  alloy,  since  it  was  the  result  of 
regression  analysis  performed  on  a  set  of  crack  growth  data  for  this  material.  According  to  reference  [7], 
equation  (4.1)  may  also  be  applied  to  7075-T6  aluminium.  In  order  to  make  this  expression  more 
generally  applicable,  and  also  to  make  use  of  the  McAir  experience  on  this  subject,  it  was  decided  to 
incorporate  correction  factor  Xj^  of  equation  (2.1a)  into  the  ARL  model  as  well.  This  factor  is  to  bring 
into  account  the  behaviour  of  the  material  under  cyclic  loading.  It  is  governed  by  the  ratio  fr  between 
the  cyclic  and  monotonic  yield  stress.  According  to  reference  [12],  this  ratio  is  approximately  equal  to 
1.12  for  both  the  2024-T3  and  the  7075-T6  aluminium  alloy.  To  obtain  a  value  of  one  for  correction 
factor  x^r  when  applied  to  these  two  materials,  this  value  of  1.12  was  used  to  normalize  the  ratio  fr.  The 
expression  for  the  stress  ratio  correction  curve  as  used  for  the  ARL  model  then  became: 

where 

X,  =  0.45  +  0.22R  +  0.21  R^  +  0.12R3 
Xf,  =  1  -  (0.756  -  0.912fr'  +  0.156fr’2)(l  -  R) 
fr'  =  fr/1.12 

The  expression  for  x^  was  simply  derived  from  equation  (4,1). 

It  may  seem  rather  strange  that  parameter  x,  has  been  included  in  equation  (4.2),  while  parameter 
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Xjj  of  equation  (1.1)  is  not  employed.  In  section  2.3.2  it  was  shown  that  the  latter  parameter  may  give 
rise  to  large  differences  in  crack  growth  rate  when  used  to  transform  thin  specimen  constant  amplitude 
data  to  thick  specimen  data,  and  vice  versa.  It  is  not  clear  whether  this  is  realistic  or  not.  Publications  on 
the  effect  of  thickness  on  crack  growth  rate  are  sometimes  contradictory,  even  when  due  to  the  same 
author  [13,19].  This  is  why  parameter  was  not  included  in  the  ARL  model,  though  a  deeper 
investigation  may  be  worthwhile. 

It  must  be  noted  here  that  the  growth  rate  curve  for  the  A7-U4SG-T651  aluminium  alloy,  described 
in  section  3.2,  was  not  corrected  for  the  cyclic  yield  stress.  The  only  reason  for  this  is,  that  the  growth 
rate  curve  was  derived  before  much  was  known  about  the  McAir  model. 

4.1.2  Variable  Amplitude  Loading 

As  for  variable  amplitude  loading,  the  ARL  model  differs  from  the  McAir  model  on  two  important 
points.  The  first  point  concerns  a  peculiarity  that  is  inherent  to  the  way  in  which  load  blocks  arc  handled 
by  the  McAir  model.  The  second  point  concerns  the  minimum  effective  stress  intensity  in  a  load  cycle  in 
relation  to  the  ovcrioad/underload  combination  that  occurred  prior  to  that  particular  load  cycle. 

The  first  point  is  best  explained  from  figure  4.1.  This  figure  shows  a  lo-hi  stress  sequence  that 
consists  of  two  load  blocks.  The  dashed  line  indicates  the  minimum  effective  stress  intensity  in  a  cycle. 
Since  the  McAir  model  operates  on  blocks  rather  than  on  individual  cycles,  the  minimum  effective  stress 
intensity  will  be  constant  within  a  block  and  the  effective  stress  intensity  range  AK^ff  of  the  first  cycle  of 
the  block  will  be  as  indicated  in  the  upper  part  of  figure  4.1.  This  behaviour  is  not  very  realistic, 
however,  for  the  material  will  not  'feel'  that  a  new  block  of  cycles  has  started  until  point  A  is  reached. 
The  situation  that  is  depicted  in  the  lower  part  of  figure  4.1  seems  more  realistic,  and  was  therefore 
adopted  for  the  ARL  model.  For  this  purpose  the  minimum  effective  stress  intensity  and  also 

the  crack  increment  Aa  for  that  matter,  is  calculated  cyclc-by-cycle  in  the  ARL  model.  The  effective 
su'css  intensity  range  AK^ff  then  is  calculated  as  the  maximum  suess  intensity  in  the  current  cycle 
minus  the  minimum  effective  stress  intensity  in  the  previous  cycle. 

This  procedure  may  seem  rather  trivial  if  very  structured  sequences  arc  considered  with  many  load 
cycles  per  block.  It  docs  not  really  matter  then  how  accurately  the  effective  stress  intensity  range  in  the 
first  cycle  of  each  block  is  calculated.  This  is  not  true  if  crack  growth  is  mainly  caused  by  overloads, 
however,  since  they  arc  most  likely  grouped  in  .small  blocks.  In  that  case  inaccurate  calculation  of  the 
effective  stress  intensity  in  the  first  cycle  may  well  affect  the  overall  result.  In  that  respect  table  6  is  of 
much  interest.  This  table  shows  the  contribution  of  overloads  to  crack  growth  for  the  load  sequences  that 
were  used  in  testing  the  McAir  model.  The  number  of  overloads  as  a  percentage  of  the  total  number  of 
cycles  is  also  listed.  For  reasons  of  convenience  this  table  was  produced  using  the  ARL  model.  This  docs 
not  matter  much,  however,  since  the  definition  of  overloads  in  both  models  is  practically  identical.  Each 
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lime  two  figures  are  listed,  indicating  a  range  of  values  that  is  covered  while  the  crack  grows.  The  first 
figure  describes  the  situation  just  after  crack  growth  started,  and  the  second  figure  is  associated  with  the 
critical  crack  length.  The  first  figure  is  more  important,  since  it  describes  about  80%  of  the  crack  growifi 
life.  For  instance,  the  contribution  of  overloads  to  crack  growth  for  the  f-b-f  sequence  with  a  stress  level 
of  1 1.09  MPa/g  started  at  9%  and  went  up  slowly  to  12%  at  80%  of  the  crack  growth  life,  and  only  after 
that  it  quickly  went  up  to  25%  at  failure.  It  may  seem  rather  strange  that  the  relative  number  of 
overloads  should  vary  with  the  crack  length.  It  must  be  realized,  however,  that  a  new  overload  is  defined 
when  the  crack  has  grown  through  the  plastic  zone  due  to  the  previous  overload.  Since  the  plasic  zone 
size  is  proportional  to  and  crack  growth  approximately  to  AX'*,  it  is  not  surprising  that  the 

relative  number  of  overloads  increases  with  increasing  crack  length.  From  table  6  it  can  be  concluded 
that  crack  growth  under  the  retardation  sequence  is  dominated  by  overloads,  while  overloads  play  a 
moderately  important  role  under  the  remaining  sequences.  Thus,  the  new  approach  of  using  the 
minimum  effective  suess  intensity  of  the  previous  cycle  when  calculating  AK^ff  seems  justified. 

The  second  point  on  which  the  ARL  model  differs  from  the  McAir  model  regarding  variable 
amplitude  loading,  is  best  explained  from  figure  4.2.  TTiis  figure  is  equivalent  to  figure  2.5,  which 
applies  to  the  McAir  model.  For  a  specific  material,  the  minimum  effective  stress  intensity  in  the  current 
cycle  is  only  a  function  of  the  previous  overload/underload  combination,  as  opposed  to  the  McAir 
model,  where  the  minimum  effective  stress  intensity  is  also  a  function  of  the  current  stress  ratio.  The 
difference  in  result  when  applying  both  models  to  a  certain  load  sequence  is  also  indicated  in  figure  4.2. 
According  to  the  McAir  model,  the  minimum  effective  suess  intensity  in  the  cycles  following  an 
overload  will  go  up  with  increasing  stress  ratio,  while  the  ARL  model  will  predict  the  minimum 
effective  stress  intensity  to  remain  constant.  It  is  not  obvious  why  the  designers  of  the  McAir  model 
have  decided  to  make  the  minimum  effective  stress  intensity  a  function  of  the  stress  ratio  of  the  emrent 
cycle.  There  is  no  apparent  physical  argument  for  this,  though  the  reason  may  have  been  to  avoid  the 
maximum  effective  stress  intensity  range  AK^jf  to  become  greater  than  the  nominal  stress  intensity  range 
AK  in  a  cycle.  In  the  ARL  model  this  is  avoided  by  including  an  'if-siatcment'  in  the  computer  code  that 
sets  AKjff  equal  to  AK  whenever  the  latter  is  exceeded.  An  exception  to  this  rule  is  the  first  cycle  of  a 
block,  as  is  explained  in  figure  4.3.  It  may  be  clear  from  comparison  of  figures  2.5  and  4.2  that  the  effect 
of  underloads  is  more  emphasised  by  the  ARL  model  than  by  the  McAir  model. 

Overloads  and  underloads  in  the  ARL  model  arc  updated  in  the  same  way  as  described  in  section 
2.3.1  for  the  McAir  model.  Two  minor  differences  arc  that  (1)  the  update  criterion  check  is  given  by  the 
dashed  line  in  figure  2.6,  and  (2)  the  size  of  the  plastic  zone  due  to  an  overload  is  given  by: 

w,,,  =  (P/7t)(K^„,„|/Oy,ff)2  (4.3) 


where 
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P  =  0.333  +  0.5a  +  0.167a2 

In  this  way  the  plastic  zone  under  plane  strain  is  three  times  larger  than  under  plane  stress.  According  to 
reference  [8],  from  which  equation  (4.3)  was  derived  (except  for  Oyeff.  which  was  obtained  from  [6]), 
this  is  more  common  than  the  factor  of  two  that  is  used  in  the  McAir  model.  The  equation  for  p  simply 
evolved  from  the  desire  that  a  value  for  a  of  zero  corresponds  to  plane  stress,  a  value  of  one  to  plane 
strain,  and  a  value  of  minus  one  to  minimum  retardation  (P=0),  as  is  also  the  case  for  the  McAir  model. 
It  should  be  realized  that  factor  a  (and  p  for  that  matter)  is  only  used  for  the  calculation  of  the  plastic 
zone  size,  while  the  McAir  model  also  applies  this  factor  in  the  determination  of  as  can  be 

appreciated  from  equation  (2.1). 

Unlike  the  McAir  model,  the  ARL  model  does  not  incorporate  a  maximum  effective  stress  intensity 
factor.  This  implies  that  sustained  crack  arrest  after  a  high  overload  can  occur,  as  far  as  the  ARL  model 
is  concerned.  Like  the  McAir  model,  the  ARL  model  does  not  distinguish  between  single  and  multiple 
overloads. 

4.2  Predictions 

Crack  growth  was  predicted  using  program  ARLECC  and  the  four  sequences  described  in  section 
3.1,  for  the  two  stress  scaling  levels  of  11.09  and  17.00  MPa/g.  Again  the  material  properties  of  tables  1 
and  2  were  used.  Given  the  rather  large  specimen  thickness  a  value  of  zero  was  used  for  input  parameter 
a. 

Predicted  crack  growths  arc  shown  in  figure  4.4  for  the  four  sequences  and  the  two  stress  levels 
examined.  In  figure  4.5  the  predictions  from  the  ARL  model  arc  compared  against  the  experimental 
crack  growth  data.  This  figure  also  shows  the  predictions  for  the  unretarded  or  no-load-intcraction  case. 
These  latter  predictions  were  made  using  the  associated  option  of  program  ARLECC.  This  option 
involved  the  use  of  cquau’on  (4.2)  to  account  for  the  stress  ratio  effect  in  each  cycle.  The  predicted  crack 
growth  lives  arc  summarized  in  tables  7  and  8,  which  arc  equivalent  to  tables  3  and  4  for  the  McAir 
model.  In  principle,  the  corresponding  values  of  the  no-load-interaction  predictions  in  tables  3  and  7 
should  be  equal.  Since  the  ARL  model  calculates  crack  growth  cycle-by-cycle,  and  because  this  model 
incorporates  a  different  stress  ratio  correction  curve,  the  results  differ  slightly  however. 

The  crack  growth  lives  were  predicted  quite  well  for  both  sacss  levels.  The  maximum  difference 
between  predicted  and  experimental  lives  is  1.28:1  for  the  random  sequence  at  a  sucss  level  of  11.09 
MPa/g.  and  1.39:1  for  the  acceleration  sequence  at  17.00  MPa/g.  Apparently  the  theoretical  value  of 
zero  for  input  parameter  a  was  very  appropriate,  hence  no  sensitivity  analysis  of  predicted  crack  growth 
to  a  was  undertaken.  The  most  important  observation,  however,  is  that  the  experimental  trend  of  the 
structured  sequences  giving  ll?c  longest  crack  growth  lives  was  predicted  correctly.  As  such,  the  ARL 


19 


model  is  superior  to  other  models  examined  so  far. 

To  check  whether  the  ARL  model  predicts  the  experimental  trend  of  thin  specimens  showing  a 
longer  delay  in  crack  growth  after  a  single  overload  than  thick  specimens,  computer  program  ARLECC 
was  used  to  analyze  the  conriguradon  that  is  described  in  section  3.4.  Again  a  single  overload  was 
considered,  followed  by  constant  amplitude  R=0  loading  with  a  maximum  cyclic  stress  of  50  MPa. 
Only  two  values  for  a  were  considered,  namely  zero  and  one. 

The  results  of  the  calculations  are  presented  in  figure  4.6.  The  corresponding  results  from  the 
McAir  model  are  also  included  in  this  figure.  Again  the  experimental  trend  is  predicted  correctly.  The 
delay  period  Nj  after  an  overload  is  approximately  three  times  longer  for  thin  specimens  (a  =  1)  than  for 
thick  ^lecimens  (a  =  0).  This  value  is  directly  related  to  the  difference  in  plastic  zone  sizes,  since  the 
stress  ratio  correction  curve  does  not  change  with  factor  a.  Unlike  the  McAir  model,  the  ARL  model 
predicts  sustained  crack  arrest  after  very  high  overloads.  The  shut-off  ratio  /<!„,,  at  which  this 
occurs  can  be  derived  from  equation  (4.2)  and  is  given  by; 

X,  Xf, 

Equation  (4.4)  is  independent  of  factor  a.  For  the  current  example  this  equation  yields  a  shut-off  ratio  of 
2.54.  According  to  figure  4.6,  however,  overload  ratios  greater  than  two  already  give  very  long  delay 
periods.  It  is  interesting  to  note  that  the  shut-off  ratio  in  figure  3.6  is  equal  to  approximately  2.7.  This 
supports  the  ARL  model,  although  a  different  material  under  slightly  different  conditions  was  considered 
in  thb  figure. 

In  the  McAir  model  the  stress  ratio  correction  varies  with  the  specimen  thickness,  whereas  this 
curve  is  independent  of  a  in  the  ARL  model.  Together  with  the  difference  in  calculation  of  the  plastic 
zone  size,  compare  equations  (2.2)  and  (4.3),  this  explains  the  difference  between  the  McAir  curves  and 
the  ARL  curves  in  figure  4.6.  The  'straight'  part  of  the  McAir  curves  was  already  explained  in  section 
3.4. 


4J  DisetiMioa 

For  both  stress  levels  examined,  iltc  predictions  from  the  ARL  model  were  quite  good.  Apparently 
the  theoretical  value  of  zero  for  input  parameter  a  was  quite  appropriate.  The  cxperimcnml  trend  of  the 
structured  sequences  giving  longer  crock  growth  lives  was  predicted  correctly.  The  effect  of  .specimen 
ihickncsf  on  spectrum  c.rock  growth,  os  dc,scribed  in  the  literature,  was  also  predicted  correctly. 

Comptrisc/n  of  ubics  3  and  7  allows  that  the  rc.sulLs  from  the  McAir  model  and  the  ARL  model 
differ  mainly  for  the  more-random  sequences,  tlte  ARl,  results  being  only  about  58  %  of  the 
corresponding  McAir  reiulu.  Tills  Is  mainly  caused  by  Ihc  more  pronounced  effect  of  underloads  on 
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crack  growth  in  the  ARL  model.  It  is  true  that  the  ARL  model  differs  on  more  points  from  the  McAir 
model,  but  it  can  be  made  plausible  that  these  points  are  not  very  important  in  explaining  the  difference 
in  results  from  the  two  models.  Firstly,  the  ARL  model  uses  the  minimum  effective  stress  intensity  of 
the  previous  cycle  to  calculate  the  effective  stress  intensity  range  in  the  current  cycle,  as  shown  in  figure 
4.1.  It  was  argued  in  section  4.1.2  that  this  is  only  important  if  crack  growth  is  mainly  determined  by 
overloads  that  occur  in  small  blocks  of  only  a  few  cycles.  It  was  shown  that  this  was  not  the  case  for  the 
random  and  the  f-b-f  sequence.  For  these  sequences,  overloads  only  determined  about  10  to  20  %  of 
crack  growth,  depending  on  the  stress  level.  This  is  not  enough  to  explain  the  large  differences  between 
the  results  from  the  tvro  models.  In  fact,  only  crack  growth  under  the  retardation  sequence  was  mainly 
determined  by  overloads,  see  table  6,  but  these  overloads  were  mainly  grouped  in  large  blocks,  as  can  be 
appreciated  from  figure  3.1b. 

Secondly,  the  ARL  model  employs  a  different  stress  ratio  correction  curve,  together  with  a 
different  equation  for  the  plastic  zone  size.  From  figure  4.6  it  follows  that  for  a  thick  specimen  the  ARL 
model  predicts  more  crack  growth  delay  after  a  single  overload  than  the  McAir  model.  This  indicates 
that  the  use  of  the  new  stress  ratio  correction  curve,  together  with  the  different  equation  for  plastic  zone 
size,  gives  rise  to  more  crack  growth  retardation  in  thick  specimens  under  variable  amplitude  loading 
instead  of  less,  as  suggested  by  the  difference  between  tables  3  and  7.  This  leaves  the  difference  in 
modelling  of  the  underloads  as  the  only  explanation  why  the  ARL  model  gave  better  results  than  the 
McAir  model.  It  was  already  shown  in  table  3  that  the  assumption  of  the  crack  being  stationary  during  a 
load  block  had  little  effect  on  the  accuracy  of  the  results. 


5.  CONCLUSIONS 

The  following  conclusions  may  be  drawn  from  the  previous  sections,  which  tested  the  McAir 
model  against  experimental  crack  growth  that  covered  several  sequences  and  two  s'rcss  scaling  levels.  It 
must  be  realized  that  these  conclusions  pertain  to  a  particular  version  of  the  McAir  model,  dating  back  to 
1984,  Later  versions,  which  unfortunately  were  not  available  to  ARL  at  the  time  this  report  was  written, 
may  have  been  drastically  revised. 

1 .  The  McAir  model  predicted  crack  growth  under  the  more-structured  sequences  quite  well. 

2.  The  morc-rand(Kn  sequences  were  handled  less  well  by  the  McAir  model.  Very  uncon-servative 
ftwiu  were  obtained.  Even  a  value  as  low  as  -  I  for  input  parameter  a  could  not  remedy  this  for  the 
lo»cf  stress  level.  This  is  caused  by  poor  modelling  of  the  effect  of  underloads  in  the  McAir  model. 

1  The  wend  in  cxpcnmcntal  crack  growth  of  Uic  morc-.struciurcd  sequences  giving  longer  lives  than 


the  more-random  sequences  was  not  predicted  correctly.  In  this  respect,  the  McAir  model  is  not  belter 
than  simpler  models  like  the  Wheeler  model  and  the  Modified  Willenborg  model. 

4.  Input  parameter  a  is  the  only  parameter  that  can  be  used  to  tune  the  McAir  model.  This  parameter 
was  shown  to  represent  the  specimen  thickness  and,  as  such,  has  some  physical  significance.  Higher 
stress  levels  tend  to  require  a  h’^ther  value  for  this  parameter,  however. 

5.  A  proposed  ARL  model,  based  somewhat  on  the  McAir  model,  predicted  crack  growth  under  the 
four  load  sequences  and  the  two  stress  levels  quite  well.  It  is  the  only  model  examined  so  far  on  these 
data  that  correctly  predicted  the  experimental  trend  of  the  more-structured  sequences  giving  the  longest 
crack  growth  lives.  For  the  basic  option  of  a  through  crack  in  a  CCT-spccimen,  the  proposed  ARL  model 
is  superior  to  the  McAir  model. 

6.  The  experimental  trend  of  thin  specimens  showing  longer  crack  growth  delay  after  a  high  overload 
than  thick  specimens  was  predicted  correctly  by  both  models.  Unlike  the  ARL  model,  the  McAir  model 
predicts  that  no  sustained  crack  arrest  will  occur  after  a  very  high  overload. 

7.  Neither  model  accounts  for  the  effect  of  multiple  overloads.  The  same  applies  to  the  phenomenon 
of  delayed  retardation  after  an  overload. 
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AK 

(MPaVm] 

daMN 

Im/cycio) 

I.4S4SS 

1.000000-10-‘' 

1.71190 

4.000000-10-‘' 

120198 

l.993473'10-'° 

3.09742 

910141610-*® 

26.9774 

2.864725-ia‘ 

33.0770 

8310194-10^ 

37.0996 

2.461269-10-S 

39.9125 

7389003- 10-5 

41.9771 

2338263-10-< 

4410S8 

1.193590  10-J 

45.7876 

6.425438- 10-3 

46.9660 

3.481965-10-2 

47.8767 

1.896I3310-' 

48.6000 

1.0 

Table  1 .  Crack  growth  curve  for  A7>U4SG>T65 1  alumimium  alloy,  stress  ratio  R  =  0. 


Young's  modulus 

69000  MPa 

Monoionic  yield  stress 

457  MPa 

Cyclic  yield  strc.ss 

411  MPa 

Prxiurc  toughness'^ 

48.6  MPaVm 

tMuccd  from  (Miguo  uiu.  A  tUlic  teil  yioMed  72.4  MPaVm. 


Table  2.  Static  properties  of  A7-LMSO-T651  alumimium  alloy. 


stress  scale 

a 

calculated  life  [programs] 

[MPa/g] 

acceleration 

retardation 

f-b-f 

random 

11.09 

0.0 

48.0 

47.0 

58.3 

57.7 

-0.2 

43.2 

43.0 

52.7 

52.3 

-1.0 

26.5 

26.9 

30.8 

31.2 

unietarded 

21.0 

21.7*) 

23.4 

23.4 

17.00 

0.0 

9.0 

8.0 

■9 

11.1 

-0.2 

8.0 

8.0 

10.0 

-1.0 

5.0 

5.9 

6.2 

unretarded 

4.0 

4.8 

4.6 

21.0  if  cakultted  cyde-by-cyck 


Table  3.  Crack  growtli  lives  as  predicted  with  the  McAir  model  for  various  values  of  a. 


stress  scale 

a 

relative  life  (%  of  experimental  life) 

[MPa/g] 

acceleration 

retardation 

f-b-f 

random 

11.09 

0.0 

119 

121 

214 

226 

-0.2 

107 

111 

194 

205 

-l.O 

65 

70 

113 

122 

unretarded 

52 

56 

86 

92 

17.00 

0.0 

81 

84 

144 

151 

-0.2 

72 

84 

130 

136 

-1.0 

45 

62 

78 

84 

unretarded 

36 

51 

59 

62 

Table  4.  Predicted  crack  growth  lives  as  a  percentage  of  the  experimental 
crack  growth  life  (McAir  model). 


stress  scale 

load  sequence 

[MPa/g] 

acceleration 

retardation 

f-b-f 

random 

11.09 

Un  (52) 

Un(56) 

Un(86) 

Un(92) 

Wh(62) 

Wh(59) 

Wh(lll) 

Wh  (122) 

MW  (71) 

MW  (69) 

MW  (144) 

MW  (154) 

Me  (107) 

Wi(98) 

Me  (194) 

Me  (205) 

Wi  (108) 

Me  (III) 

Wi(235) 

Wi  (268) 

17.00 

Un(36) 

Un(51) 

Un  (59) 

Un  (62) 

MW  (63) 

Wh(62) 

Wh(125) 

Wh  (132) 

Wh(64) 

MW  (70) 

Me  (130) 

Me  (136) 

Me  (72) 

Me  (84) 

MW  (138) 

MW  (149) 

Wi(76) 

Wi  (88) 

Wi  (169) 

Wi  (180) 

Un  no-loid-inteniction 

Wh  Wheeler  model 

Wi  Willenborg  model 

MW  Modified  Willenborg  model 

Me  McAir  model  (a  =  -0.2) 


Table  5.  Comparison  of  predicted  crack  growth  lives.  Figures  indicate  percent  of  experimental 
crack  growth  life.  Data  obtained  from  reference  [1], 


stress  scale 

load  sequence 

[MPa/g] 

(%] 

acceleration 

retardation 

f-b-f 

random 

11.09 

^ol^toul 

0.3 -0.6 

0.3  -  0.5 

0.3  -  0.5 

^o|/^^oul 

10-35 

9-25 

9-23 

17.00 

^ol/^ioul 

0.7  -  1.1 

3.5  -  4.8 

0.6  -  0.8 

0.5  -  0.9 

'^ol/'^ul 

19-39 

74-87 

14-31 

16-39 

N 

ol  toul  number  of  overloids  already  experienced 

V 

u>lal  total  number  of  cycles  already  experienced 
^ol  total  crack  growth  due  to  overloads 

**total  total  crack  growth 

Table  6.  Contribution  of  overloads  to  crack  growth  (calculated 
using  ARL  model). 


stress  scale 

a 

calculated  life  [programs] 

[MPa/g] 

acceleration 

retardation 

f-b-f 

random 

11.09 

0.0 

41.0 

34.2 

32J 

unretarded 

21.0 

mam 

23.6 

23.5 

17.00 

0,0 

8.0 

8.0 

6.8 

6.5 

unretarded 

4.0 

4.0 

4.7 

4.7 

Table  7.  Crack  growth  lives  as  predicted  with  the  ARL  model. 


stress  scale 

a 

relative  life  (%  of  experimental  life) 

[MPa/g] 

acceleration 

retardation 

f-b-f 

random 

11.09 

0.0 

101 

101 

126 

128 

unretarded 

52 

54 

87 

92 

17.00 

0.0 

72 

84 

87 

87 

unretarded 

36 

42 

60 

63 

Table  8.  Predicted  crack  growth  lives  as  a  percentage  of  the  experimental  crack 
growth  life.  The  final  life  was  used  for  comparison  (ARL  model). 
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Figu*'''  2.2  McAir  stress  ratio  correction  curves  as  a  function  of  retardation  parameter  a,  for  a 
ratio  fr  between  cyclic  and  monotonic  yield  stress  of  1.0. 
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Figure  2.3  Variation  of  the  McAir  stress  ratio  correction  curves  with  ratio  fr  between  cyclic 

and  monotonic  yield  stress. 
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Figure  2.5  (continued). 


Figure  2.7 


Simple  now  diagram  of  die  Me  Air  model,  covering  the  case  of  a  through 
crack  in  a  CCT-specimen. 
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Figure  3.3  Crack  growth  curve  for  A7-U4SG-T65 1  aluminium  alloy,  stress  ratio  R  =  0. 
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Figure  3.4  Crack  growth  as  calculated  with  the  McAir  model  for  the  various  load  sequences 

and  stress  levels  (a  -  -0.2). 
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Figure  3.5a 


Comparison  of  the  predictions  from  the  McAir  model  (a  =  -0.2)  against  the 
experimental  results  for  a  stress  level  of  1 1 .09  MPa/g. 
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Figure  3.5b  Comparison  of  the  predictions  from  the  McAir  model  (a  -  -0.2)  against  the 
cxpcrimcnuil  results  for  a  sircss  level  of  17.00  MPa/g. 
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Figure  3.6  Delay  period  Nj  as  a  function  of  the  overload  level  in  7075-T6  aluminium, 

found  experimentally  [15]. 
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Figure  4.3  Calculation  of  the  minimum  effective  stress  intensity  under  variable 

amplitude  loading  (ARL  model). 


Figure  4.3 


The  effective  stress  intensity  range  AKj,ff  of  the  first  cycle  of  a  block  may  be  greater 
than  the  nominal  stress  intensity  range  AK  (ARL  model). 


Figure  4.4  Crack  growth  as  calculated  with  the  AI 
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Figure  4.5a  Comparison  of  the  predictions  from  the  ARL  model  against  tlic  experimental 
results  for  a  stress  level  of  1 1.09  MPa/g. 
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Figure  4,5b 


Comparison  of  the  predictions  from  the  ARL  model  against  the  experimental 
results  for  a  stress  level  of  17.00  MPa/g. 
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APPENDIX  A  COMPUTER  PROGRAM  CNT8SHR 

This  appendix  describes  iho  computer  program  CNT8SHR.  This  program  incorporates  the  McAir 
model  as  described  in  this  repott.  It  is  a  stripped  version  of  computer  program  CNTKM8  that  was 
obtained  from  McDonnell  Douglas.  According  to  a  comment  in  the  source  listing  this  latter  program 
dates  back  to  November  1984.  It  is  the  latest  release  available  to  ARL.  Apart  from  the  user's  manual  the 
code  of  CNTKM8  was  the  only  source  that  provided  detailed  information  on  the  McAir  model.  To  be 
able  to  extract  the  necessary  information  it  proved  useful  to  rewrite  the  original  source  code,  since  this 
code  had  been  written  in  a  rather  unstructured  way.  Furthermore  the  original  program  covered  many 
options  that  are  not  essential  for  the  understanding  of  the  basic  McAir  model.  This  finally  resulted  in 
computer  program  CNT8SHR.  This  program  only  covers  the  case  of  a  through  crack  in  a  CCT 
specimen.  Care  has  been  taken  that  the  modified  version  and  the  original  program  are  fully  compatible, 
so  that  the  same  input-files  can  be  used.  It  has  been  cheeked  that  the  same  results  arc  obtained  for  both 
programs. 

The  computer  program  CNT8SHR  has  been  written  in  Fortran  77.  It  can  be  run  on  a  personal 
computer  with  a  mathematical  co-processor.  Small  jobs  with  simple  load  spectra  will  take  up  to  a  few 
minutes,  while  in  the  case  of  a  flight-by-flight  spectrum  a  job  may  run  for  more  than  an  hour.  The  input 
should  be  in  Sl-units  without  prefixes.  The  flow  diagram  of  CNT8SHR  is  given  in  figure  A.l. 

A.l  Interpolation  between  data  points  on  crack  growth  curve 

Both  the  original  and  the  modified  version  of  the  computer  program  expect  the  material's  R=0 
crack  growth  curve  to  be  input  in  ubuiated  form.  The  way  the  two  programs  interpolate  between  the 
data  points  differs,  however.  To  facilitate  interpolation  Uic  original  program  CNTKM8  tries  to  linearize 
the  da/dN-AK  curve  by  applying  the  following  transformation  to  the  input  data: 

X  =  logfAKjff) 

K.K 

y  =  log(  -*'j;  -da/dN) 

wlicrc  is  the  matcrial'.s  fracture  toughness  and  the  effective  stress  intensity  range  as  calculated 
with  equation  (2,1).  The  transformation  of  the  da/dN-dala  is  obviously  based  on  Forman's  equation  for 
constant  amplitude  fatigue  crack  growth.  This  u-ansformaiion  may  .seriously  affect  the  accuracy  of 
interpolation,  however.  This  is  especially  Uic  ca.sc  when  the  su-aighl  'Paris-parl'  of  the  crack  growth 
curve  IS  represented  by  the  end-points  only  (which  is  enough  anyway).  This  can  Ik  illustrated  with  the 


A,2 


following  example.  Consider  ihc  ncliiious  though  realistic  R  =  0  crack  growtli  curve  of  which  tlic 
straight  part  is  given  by: 

logfdaAdN)  »  -31.75  +  3.5 log(AK)  (A.2) 

with  daAdN  in  m/cyclc  and  AK  in  PaVm.  Assuming  that  =  50  MPaVm  and  that  AK^f  is  equal  to  50% 
of  AK.  this  curve  could  be  tabulated  as  follows: 


log(AK) 

log(da/dN) 

>< 

y 

6.50 

-9.00 

6.20 

-9.028 

7.50 

-5.50 

7.20 

-5.935 

Now  consider  a  load  cycle  that  results  in  an  effective  stress  intensity  range  of  5  MPaVm,  while  the 
nominal  stress  intensity  range  is  10  MPaVm  and  R  =  0.  With  equation  (A.l)  it  is  found  then  that 
X  =  6.70  and  linear  interpolation  in  the  table  above  yields  y  =  -7.482.  Using  equation  (A.l)  again  it  is 
found  that  the  crack  growth  rate  da/dN  is  equal  to  4. 13' 10  *  m/cycle.  Using  equation  (A.2),  however,  it 
is  found  that  da/dN  is  equal  to  5.62- 10  *  m/cycle.  It  is  clear  that  the  latter  result  is  the  right  one. 

The  difference  between  the  two  da/dN-values  is  36%.  It  is  thus  concluded  that  the  Forman-like 
transformation  may  suffice  for  very  high  AK  values,  with  near  the  fracture  toughness,  but  not  for 
moderate  values  at  all.  It  was  therefore  decided  not  to  u.se  this  transformation  in  the  modified  version 
CNT8SHR  of  the  original  program.  It  is  left  to  the  user  to  discretize  the  crack  growth  curve  in  such  a 
way  that  the  extremes  of  the  curve  are  sufficiently  accounted  for.  An  example  is  given  in  figure  3.3. 

It  should  be  noted  here  that  the  Forman-like  transformation  is  still  applied  by  the  original  program 
CNTKM8.  If  this  program  is  used  this  fact  should  be  taken  into  account.  It  may  be  necessary  then  to 
increase  the  number  of  data  points  that  are  used  to  represent  the  'Paris-part'  of  the  crack  growth  curve. 


A.2  Example  problems 

The  input  consists  of  three  data  files  and  a  straightforward  interactive  part.  The  input  should  be  in 
Sl-units  without  prefixes,  though  this  is  only  important  in  the  ease  of  artificially  cracked  specimens 
without  a  load  history.  In  that  ease  no  residual  plastic  zones  arc  present  around  the  crack  tip  just  before 
cyclic  loading  starts,  and  the  minimum  effective  stress  intensity  in  the  first  cycle  is  equal  to  zero.  The 
program,  however,  always  assumes  that  a  pla.stic  zone  is  present  at  the  time  of  the  first  cycle,  based  on 
.some  overload  that  occurred  before  analysis  starts.  Normally  this  overload  is  assumed  to  be  equal  to  the 
highest  load  in  the  input  sequence,  but  for  machined  slots  this  overload  would  be  zero.  To  avoid 
singul.arity  problems,  however,  this  overload  is  set  to  the  arbitrary  small  value  of  1000  Pa.  This  is 


automaiicaliy  done  by  the  program.  If  imperial  units  (Icsi)  or  Sl-uniis  with  prefixes  (MPa)  arc  used,  a 
value  of  1000  is  very  high,  and  this  is  why  the  input  should  be  in  Sl-uniis  without  prefixes  in  ease  of 
machined  slots.  If  fffc-cycled  specimens  are  considered  any  consistent  set  of  units  can  be  used.  It  should 
be  noted  here  that  the  original  program  assumes  imperial  units,  although  the  same  argument  as  above 
holds. 

The  three  data  files  to  be  created  arc:  The  analysis  file,  the  material  file  and  the  spectrum  file. 
These  files  will  be  explained  from  an  example  problem.  The  thick  CCT  specimen  that  is  considered  is 
made  from  A7-U4SG-T651  aluminium  alloy.  The  properties  of  this  material  are  given  in  tables  1  and  2 
of  this  report  The  width  equals  100  mm,  and  the  initial  machined  crack  is  prc-cycled  to  a  total  length  of 
40  mm.  A  blocked  load  program  is  applied  with  a  maximum  gross-area  stress,  or  design  limit  stress 
DLS,  of  60  MPa.  The  load  program  (=  spectrum  file)  consists  of  5  blocks  of  load  cycles,  and  covers  100 
flight  hours.  Two  cases  will  be  considered,  the  first  example  being  straightforward  while  in  the  second 
example  a  high  overload  of  120  MPa  is  applied  before  starting  the  fatigue  program. 

A.2.1  Analysis  file 

The  analysis  file  that  covers  this  example  is  given  in  table  A.!.  It  is  named  example. anl.  This  file 
conuins  the  information  that  is  not  covered  by  the  other  files  or  by  the  interactive  input. 

The  first  line  is  a  commentary  line.  It  can  contain  up  to  72  characters  and  should  not  be  empty. 
The  second  line  contains  an  integer  dummy  variable.  This  variable  is  only  used  by  the  original  program 
CNTKM8.  The  third  line  contains  the  specimen  width,  or  twice  the  shortest  distance  from  the  plate  edge 
to  the  crack  centre  in  case  of  asymmetric  cracks.  Lines  4  and  5  contain  integer  dummy  variables  that  are 
only  used  by  the  original  program.  The  sixth  line  gives  the  material  file  name.  The  seventh  line  gives 
the  maximum  number  of  load  programs  to  be  analyzed.  This  is  the  cut-off  for  very  slow  crack  growth. 
The  maximum  number  is  400.  This  number  is  determined  by  the  dimensions  of  the  arrays  used  by  the 
program.  The  last  line  contains  the  number  of  hours  between  damage  print-out  It  is  advisable  to  input  a 
number  that  is  an  integer  multiple  of  the  number  of  hours  that  is  covered  by  the  load  program. 

A.2.2  Material  file 

The  material  file  designates  the  da/dN-AK  data  for  use  in  the  analysis.  The  program  transforms  the 
data  to  a  da/dN-AK^jf  curve,  assuming  that  the  input  curve  is  for  a  stress  ratio  of  zero.  The  material  file 
used  in  this  example  is  given  in  table  A.2,  and  is  named  dadn.dat. 

The  first  line  of  the  materia!  file  is  a  commentary  line.  It  may  contain  up  to  72  characters  and 
should  not  be  empty.  The  second  line  contains  Young’s  modulus,  Poisson's  ratio,  the  monotonic  yield 
stress  and  the  cyclic  yield  stress.  It  should  be  noted  here  that  Poisson's  ratio  actually  is  a  dummy 
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parameter  for  CNT8SHR.  It  is  only  used  by  the  original  program.  The  third  line  contains  the  fracture 
toughness  and  the  integer  number  of  points  in  the  da/dN-AK  uble.  This  number  should  not  be  higher 
than  50.  To  avoid  singularity  problems,  the  fracture  toughness  must  be  slightly  higher  than  the  highest 
AK  that  occurs  in  the  table.  The  remaining  lines  contain  the  actual  da/dN-AK  table. 

A.2.3  Spectrum  file 

Spectra  (strictly:  sequences)  may  be  of  any  length  and  form:  Random,  blocked  and  flight-by-flighi. 
The  stresses  are  input  as  a  proportion  of  the  design  limit  stress,  which  serves  as  a  kind  of  stress  scaling 
factor.  The  design  limit  stress  is  input  manually.  In  that  way  it  is  easy  to  perform  a  parametric  study 
involving  the  stress  level.  For  this  example  two  stress  spectra  are  created.  One  for  the  'regular'  case, 
named  exempli. seq,  and  one  for  the  overload  case,  named  exampl2.seq.  They  are  given  in  tables  A.3a 
and  A.3b.  The  two  files  are  identical,  except  for  the  third  line. 

Again,  the  first  line  of  the  spectrum  file  is  a  commentary  line.  It  may  contain  up  to  72  characters 
and  should  not  be  empty.  The  second  line  contains  the  (integer)  number  of  blocks  nlevs  of  load  cycles  in 
the  program,  and  the  number  of  flight  hours  that  is  represented  by  the  program.  For  random  or  cycle-by¬ 
cycle  sequences  each  block  would  consist  of  one  cycle. 

The  third  line  contains  the  highest  and  lowest  stress  previously  experienced,  as  a  proportion  of  the 
design  limit  stress.  These  stresses  are  used  to  initialize  Usually  these  stresses  are  equal  to  the 

highest  and  lowest  that  occur  in  the  sequence.  In  the  case  of  a  machined  slit  that  has  not  been  pre-cycled 
these  stresses  must  be  set  to  zero.  This  tells  the  program  that  no  residual  plastic  zones  are  present  around 
the  crack  tip  yet.  It  is  also  possible  to  input  a  higher  suess  than  present  in  the  program.  Effectively  this 
means  the  application  of  a  high  overload  just  before  cycling  starts.  Comparison  of  tables  A.3a  and  A.3b 
shows  that  in  the  second  example  an  overload  equal  to  twice  the  design  limit  stress  is  applied,  as  was 
required. 

The  next  nlevs  lines  contain  the  actual  load  specuiim.  Each  line  contains  a  block  of  cycles  in  the 
following  format:  Maximum  cyclic  stress,  minimum  cyclic  stress,  number  of  cycles  in  block.  Again,  the 
stresses  are  input  as  a  proportion  of  the  design  limit  stress.  The  crack  length  is  assumed  to  be  stationary 
within  a  block  of  cycles.' This  means  that  the  stress  intensity  is  assumed  to  be  constant  during  a  block.  It 
therefore  is  advisable  to  limit  the  number  of  cycles  within  a  block  to  avoid  inaccurate  answers  in  case  of 
highly  loaded  specimens.  It  may  be  necessary  then  to  split  a  block  up  into  two  or  more  pans.  The  last 
line  contains  an  integer  dummy  parameter, 

A.2.4  Interactive  input 


A  few  parameters  have  to  be  input  by  hand.  This  facilitates  parametric  studies  involving  these 
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variables.  The  input  is  interactive  and  is,  for  the  greater  pan,  self-explanatory.  Tables  A.4a  and  A.4b 
show  what  the  screen  looks  like  when  running  the  two  examples.  The  user's  input  is  indicated  in  italic. 
Successively  the  program  asks  for  The  name  of  the  analysis  file,  the  name  of  the  spectrum  file,  the 
value  of  retardation  parameter  a,  the  initial  half  crack  length,  the  design  limit  stress,  the  amplification 
factor,  and,  after  the  program  has  finished,  the  name  of  the  output  file. 

The  retardation  parameter  ot  theoretically  accounts  for  the  thickness  of  the  specimen.  A  value  of 
one  would  pertain  to  the  situation  of  plane  stress  around  the  crack  tip,  as  present  in  very  thin  specimens, 
and  a  value  of  zero  would  pertain  to  the  situation  of  plane  strain,  as  in  the  case  of  very  thick  specimens. 
In  practice  this  factor  a  is  mainly  used  as  a  tuning  parameter.  For  this  purpose  its  value  can  range 
between  -1  (minimum  retardation)  and  -►1  (maximum  retardation). 

The  amplification  factor  is  related  to  the  failure  of  the  specimen.  Failure  is  predicted  when  some 
reference  stress  intensity  factor  exceeds  the  fracture  toughness.  This  reference  stress  intensity  factor  is 
calculated  using  the  current  crack  size  and  the  highest  load  in  the  spectrum,  including  the  incidental 
high  overload  that  may  have  occurred  before  cycling  started  (first  entry  on  line  3  of  spectrum  file 
greater  than  one).  If  it  is  desirable  to  relate  failure  to  the  highest  stress  in  the  program,  rather  than  to  tlie 
incidental  high  overload  before  cycling  started,  the  amplification  factor  must  be  made  equal  to  the  ratio 
of  the  highest  spectrum  load  and  the  high  incidental  overload.  In  the  case  of  the  second  example  this 
requires  an  amplification  factor  of  0.5  to  be  input.  In  absence  of  a  high  initial  overload  the  amplification 
factor  is  usually  equal  to  one.  On  the  other  hand,  an  amplification  factor  greater  than  one  may  also  be 
required  in  certain  situations.  The  USAF  Mil-Specs  (MIL-A-83444)  state  that  the  load  level  to  be  used 
for  failure  computation  can  be  higher  than  that  in  the  analysis  spectrum.  This  is  done  to  account  for 
excess  usage  severity  of  individual  aircraft  and  the  inspcciibility  of  the  location  in  question  [2]. 

The  output  files  for  the  two  examples  arc  listed  in  tables  A.5a  and  A.5b. 


A,3  Listing  of  CNT8SHR 


The  variables  labelled  with  an  asterisk  are  input  by  the  user  or  read  from  a  data  file.  The  variable 
names  are  the  same  as  used  in  the  original  program.  Only  the  number  of  variables  has  been  reduced,  due 
to  removing  a  number  of  options.  The  output  has  drastically  changed  compared  with  the  original 
program. 


PROGRAM  CNT8SHR 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c  c 
c  stripped  version  of  cntkmB,  that  docs  not  use  the  forman-  c 
c  transformation  {kc-kmax) /)tc  of  the  dadn-curve.  Calculation  of  c 
c  YLDl  placed  outside  if-block.  c 
c  Marcel  Bos,  june  1991  c 
c  c 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 


i ntege r  variables: 


ITCYC 


total  nusnbor  of  cylos  that  spaclmen  already  haa  c 

experienced 

KK  counter  that  keepa  track  of  the  number  of  timea  c 

that  intermediate  reauita  are  atored  in  c 

arraya  M  and  HI  c 

*NBLXS  max.  number  of  apectrum  biocka  to  be  analyzed  c 

*NDAOH  number  of  points  in  da/dn  table  (R-0  data)  c 

’NLCVS  number  of  load  levels  in  spectrum  c 

c 

real  variables:  c 

c 

A  current  crack  length  c 

Al(j)  intermediate  craeklength,  associated  with  HI{j)  c 

•Ain  initial  crack  length  in  [ml  c 

•ALPHA  retardation  parameter  c 

•AMP  amplification  factor,  ratio  of  highest  load  in  c 

magnified  inspection  interval  to  highest  load  c 

in  spectrum  block  c 

AOL  crack  length  at  time  of  overload  c 

BETA  parameter  that  is  derived  from  ALPHA  c 

BETAl  parameter  that  is  derived  from  ALPHA  and  XNU  c 

•BHRS  number  of  hours  represented  by  spectrum  block  c 

CHECK  critericn  used  to  determine  if  KMAXOL  should  be  c 

updated  c 

CODMAX  crack  tip  opening  at  time  of  overload  c 

CODX  some  crack  tip  opening  c 

CORFAC  correction  factor  that  is  applied  to  K  to  account  c 

for  finite  width  of  specimen  c 

•CYCLES  number  of  cycles  in  block  c 

•CYLD  cyclic  yield  stress  in  N/m2  c 

DADN  actual  crack  growth  da/dn,  associated  with  DKEFF,  c 

as  found  from  interpolation  in  modified  material's  c 
da/dn  vs.  delta-K  table  c 

•DADNB(j)  da/dn  part  of  material's  da/dn  vs.  delta-K  table  c 

for  R«0  (j.le.NDADN)  in  (m/cycle|  c 

•DELKB(j)  delta-K  part  of  material's  da/dn  vs.  delta-K  table;  c 

(later  on  converted  into  some  effective  value)  c 

DELKB(NDADN)  should  be  less  than  KC;  in  (Pa  sqrt (m) )  c 
DKEFF  equal  to  KMXEFF-KMNEFF  c 

•DLS  design  limit  stress  in  [N/mZ]  c 

•E  modulus  of  elasticity  in  [N/m2)  c 

EFFKC  part  of  equation  for  KMNEFF  c 

EFFKR  if  SR-0  then  EFFKR  -  EFFKC/Kmax  c 

FR  ratio  CYLD/TYLD  c 

HI(j)  intermediate  result,  couples  the  total  number  of  c 

experienced  hours  to  AI(j)  c 

HOURS  equal  to  Icompleted  blocks'BHRS  (before  failure)  c 

or  the  time  to  failure  o 

•Kc  material  fracture  tougness  in  (Pa  sqrt (m) )  o 

KMAX  max  K  in  cycle  c 

KMAXOL  max  K  in  overload  cycle  c 

KMIN  min  K  in  cycle  c 

KMINOL  underload  c 

KMNEFF  effective  min.  K  in  cycle  c 

KMXEFF  effective  max.  K  in  cycle  c 

PHOURS  gives  the  time  at  which  next  intermediate  results  c 

should  be  printed  c 

•PRTHRS  number  of  spectrum  hours  between  damage  printout  c 

(must  be  .GE.  BHRS)  c 

RC  delta-K/RC  gives  the  effective  K  for  R=0  c 

this  parameter  mainly  depends  on  ALPHA  c 

REFK  equal  to  AMP*SMAX1*CORFAC^DSORT (PI*A)  c 

if  REFK  exceeds  Kc  then  failure  is  assumed  c 

REFKIA  used  in  final  calculation  of  A  at  failure  c 

•SIGmax  -  max  .stress  in  r:ycie  c 

•SIGmin  min  stress  in  cycle  c 

Smax  either  equal  to  SMAXl  or  to  1000.  («=  small  number  c 

if  Sl-units)  if  specified  spectrum  file  ISFID  c 

Smaxl  highest  stress  in  spectrum  (or  higher  if  specified  c 

in  ISFID)  c 

Smin  lowest  stress  in  spectrum  (or  lower),  or  0  if  c 

Smax  “  1000.  c 

SR  stress  ratio  Kmin/Kmax  c 

•Tyld  monotonic  yield  stress  in  N/m2  c 

*W  twice  the  shortest  distance  from  plate  edge  to  c 

crack  center  In  [m]  (width)  c 

WOL  plastic  zone  size  due  to  overload  c 

X2  storage  for  A  c 

XALPHA  see  equation  (2.1)  of  report  c 

XCYC  number  of  cycles  in  a  spectrum  block  c 

XFR  see  equation  (2.1)  of  report  c 

XINCR  u.sed  in  calculation  of  A  at  failure  c 

•X)JU  Poisson'sratio  c 


n  nn 


c 

XONL 

used  in  calculation  of  A  at  failure 

c 

c 

XR 

see  equation  (2.1)  of  report 

a  s)>ectruffl  blocX 

c 

c 

YCYC 

number  of  cycles  already  covered  of 

c 

c 

to  0  when  a  block  has  baen  completed) 

c 

c 

YLDl 

sort  of  effective  yield  stress 

c 

c 

c 

c 

character 

variables: 

c 

c 

c 

c 

'lAFID 

analysis  file  name  Cold*) 

c 

c 

IBIDX 

temporary  file  that  contains  stress 

spectrum  ('new') 

c 

c 

'IPRED 

prediction  file  ('new') 

c 

c 

'ISFID 

general  spectrum  file  ('old') 

c 

c 

'MATFIL 

material  file  Cold') 

c 

cccecccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

IMPLICIT  NONE 

INTEGER  DUMMY, IFALl, IFLAG, i, II. ITCYC, J, KK, M, N, NBLKS, NDADN, NLEVS 
RFJVL*a  A,AI  taOO)  , Ain, ALPHA, AMP, AOL,BETA,BETAl,BHR5, CHECK, 

♦  CODMAX, CODX, CORFAC, CYCLES, CYLO, OADN,  DADNB ( 50 ) , hal fpi , 

♦  DELKB(50) ,DKEFF,DLS,E,EFFKC,EFFKR,FR,HI (400) , HOURS, 

♦  KC,  KMAX, KKAXOL, KHIN, KMINOL, KMNEFF, KMXEFF, PHOURS, PI , 

♦  PRTHRS, RC, REFK, REFKIA, SIGmax, SIGmin, Smax, Smaxl , Smin, 

♦  SR, Tyld, W, WOL, X2, XALPHA, XCYC, XFR, XI NCR, XNU, XONE, XR, 

♦  YCYC,YLD1 
CHARACTER'S  IBIDX 

CHARACTER*30  lAFID, IPREO, ISFID, MATFIL 
CHARACTER* 72  TITLES, TITLE 

DATA  PI,halfpi  /:.  1 415926536,  1 . 5707963268/ 

ITCYC  -  0 
IFALl  -  0 
IFLAG  -  0 
XCYC  -  0 . 

YCYC  -  0 . 

KK  -  1 

INPUT  MATERIAL  AND  ANALYSIS  DATA 
write  (*,  ' (/lx, a) ' ) 

♦  'Modified  version  of  the  HcAir  crac)t  growth  program  CNTKM8.' 
write  (•,'(lx,a)')  'No  Forman-transformation  of  the  da/dn-curve' 
write  (• , ' (lx, a/I ' ) 

♦  '  Marcel  Bos,  5/3/91.' 

write  (*, '  (lx, a)  '  I 

♦  'WARNING:  In  case  of  pre-cycled  specimens  it  does  not  matter' 
write  (*, '  (lx, a)  '  I 

♦  '  whether  Sl-unlts  or  Imperial  units  are  used,  as  long' 

write  (♦, '  (lx, a)  '  I 

♦  '  the  input  is  consistent.  In  case  of  saw-cut  specimens' 

write  (*, ' (lx, a) ' ) 

•*  '  (line  3  of  spectrum  file:  0,0)  only  Sl-units  should' 

write  (♦ ,  '  (lx,  a)  '  I 

♦  '  be  used  however,  although  the  error  that  is  introduced' 

write  (*,  '  (lx,a/)  ' ) 

♦  '  when  using  imperial  units  most  liltely  is  very  small.' 

6669  write  (• , '  (/lx, a, 5 ) ' )  'Input  ANALYSIS  file:  ' 
read  (*, '  (  a  ) ' )  lAFID 

OPEN  (UNIT-2, FILE-IAFIO, ERB-6669, STATUS- ' OLD ' ) 

6666  write  ( * , ' ( /lx, a, 5 ) ' )  'Input  SPECTRUM  file:  ' 
read  (*,  '  (  a  ) ' )  ISFID 

OPEN  (UNIT-4,FILE-ISFID,ERB-6666, STATUS- ' OLD ' ) 

WRITE  (*,280) 

road  (2,  '(  a)')  TITLES 

write  (*,  '(lx,a)')  TITLES 

READ  (2,*)  dummy 
c  . . .used  to  be  icor 


write  ( * 
read  (* 
write  (* 
read  (" 
write  ( * 
read  (• 
write  {• 
read  (• 


(/lx,a,$)')  'Input  ALPHA 


(  lx, a, 5) ' ) 
*  ) 

(  1 X ,  a ,  0 )  '  ) 
*  ) 

(  1 X ,  a ,  S )  '  ) 


alpha 

AINITIAL 

aln 

DLS 

dls 

AMP  FACTOR 
amp 


IF  (OL.-,  .ie.  0.)  THEN' 


print  *,  '  DUS  should  b«  qrascer  than  0!' 
print  ’  The  proqrem  STOPS* 

■top 
END  IF 

Hid)  -  0. 

Aid)  •  AIN 
A  >0. 

READ  12,*)  N 
READ  (2,*)  dummy 
READ  (2,*)  dummy 

...used  to  be  npta2  snd  npts3,  resp. 

read  (2, '(a)*)  MATFIL 
read  (2,  *  )  NBLKS 

road  (2,  •  )  PRTHRS 

CLOSE  (2) 

OPEN  (UNIT-X.FILE-HATFIL.STATUS-'OLD*) 

read  d,*(a)*)  TITLE 

read  d,  •  )  E,  XNU,  TYLD,  CYLD 

read  d,  *  )  KC.NDADN 

IF  (Kc  .le.  0.0)  THEN 

print  *,  '  Kc  should  be  greater  than  0!* 
print  ♦,  '  The  program  STOPS' 

stop 
END  IF 

DO  :  -  l.NDADN 

READ  (1,*)  DELKBd)  ,DADNB(I) 
end  do 

CLOSE  (1) 


READ  STRESS  EXCEEDANCE  DATA... 

Ibidx  -  'ibidOO* 

OPEN  (UNIT-3 , FILE- IBIDX, STATUS- ' NEW* , FORM- ' UNFORMATTED  * ) 

read  (4,  *  (a)  ' )  TITLE 
read  (4,  •  )  NLEVS.BHRS 

read  (4,  •  )  SMAX.SMIN 

IF  (SMAX.LT. .01)  IFLAG  -  1 

SMAX  -  SMAX'DLS 
SMIN  -  SMIN'DLS 

DO  I  -  l.NLEVS 

READ  (I,*)  SIGMAX, SIGMIN,  CYCLES 
XCYC  -  XCYC+CYCLES 
SIGMAX  »  SIGMAX’DLS 
SIGMIN  -  SIGMIN*DLS 
IF  (SIGMAX. GT. SMAX)  SMAX  -  SIGMAX 
IF  (SIGMIN.lt. SMIN)  SMIN  •  SIGMIN 
WRITE  (3)  SIGMAX, SIGMIN, CYCLES 
end  do 


READ  (4,*)  dummy 
...used  to  be  num 

CLOSE  (4) 

CLOSE  (3) 

OPEN  ( UNI T-3, FILE -IBIDX, STATUS- ' OLD’ , FORM- *  UNFORMATTED* ) 
5MAX1  -  SMAX 

IF  (ifl/>g.eq.1)  then 

i ;  SI -uni t s  then .  .  . 

SMAX  -  1000. 

else  smax  should  be  ♦/-  0.001 
SKIN  -  0. 

END  IF 

COMPUTE  CON.STANT  MODEL  PARAMETERS 

BETA  «  ALPHA*  (1  . -AI.rilA) /?. 

PF.TAl  -  ALPHA*  (1  . -AI.PIIA)  '  (1  . -XNU**2) /P. 


FR  -  CYLD/TYLD 

YLDl  -  h«lfPI*SrtAX/DSQRT(2./DCOS(h«lfPI*SMAX/TYLD)-2.) 
EFFKR  -  0.46733*DS0RT(BETA) 

EFFKR  •  EFFKR* (1. +0.6* (FR-1. )-0. 156* (FR-1.>»*2) 

RC  -  1. /(I. -EFFKR) 

DO  1  -  l.NDADN 

old:  DADNB(l)  -  DLOGIO (DADNB(l) » (KC-DELKB(I) ) /KC) 
dadnb(i)  dloglO  (dadnb (i) ) 

DELKB(I)  -  DLOGIO )DELKB(I)/RC) 

IF  C.GT.l)  then 

if  (DADNB(I) .LT.DADNB(I-l))  DADNB(I)  -DADNB(I-l) 
end  1  f 
end  do 

INITIALIZE  KEFF  CURVE  BASED  ON  HIGHEST  AND  LOWEST  LOADS  IN 
THE  SPECTRUM 

A  -  AIN 

CALL  KCOR(a,w,CORFAC) 

KMAX  -  SMAX*CORFAC*DSQRT(PI*A) 

KMIN  «  SMIN*CORFAC*DSQRT<PI*A) 

SR  «  KMIN/KMAX 

XR  »  0. 46733+0. 29401*SR+0.23866*SR**2 
XALPHA  -  1.- (l.-DSQRTIBETA) ) * (1,-SR) 

XFR  -  1 (0.756-0. 912*FR+0.156*FR**2) * (1. -SR) 

EFFKC  «  XR*XALPHA*XFR*KMAX  -  0.225*KMIN 
IF  (SR.LT.O.)  EFFKC  -  EFFKR*KMAX*DEXP (0.08*SR) 

KMAXOL  -  KMAX 

WOL  »  PI/8.* (KMAX/YLD1)**2*BETA 
CODMAX  -  BETAl*KMAX**2/2./E/YLDl 
AOL  -  AIN 
PHOURS  «  PRTHRS 


BEGIN  BLOCK  BY  BLOCK  CRACK  GROWTH  ANALYSIS 

WRITE  (*,  339)  AIN 

DO  M  »  1,NBLKS 

DO  200  N  -  l.NLEVS 

READ  (UNIT-3)  SIGMAX , SIGMIN, CYCLES 
YCYC  =  YCYC+CYCLES 
ITCYC  -  ITCYC+INT (CYCLES) 

CALL  KCOR(a,w,CORFAC) 

CORFAC  -  CORFAC*DSORT(PI'A) 

KMAX  -  SIGMAX'CORFAC 
KMIN  -  SIGMIN*CORFAC 

CHECK  AT  EACH  LOAD  LEVEL  TO  SEE  IF  LOWER  LOAD  LEVEL  IS 
AN  OVERLOAD 

IF  (SIGHAX.ge. 0.0001)  then 

IF  (sigmax  .ge.  tyld)  THEN 
WRITE  (6,*)'  ' 

WRITE  (6,*)'  ••*•*  WARNING  STRESS  '.SIGMAX, 

•  IS  ABOVE  THE  YIELD' 

WRITE  (6, *) '  ' 

WRITE  (6,*) 

'  THIS  CAUSES  A  CRASH  THERE  PROGRAM  WILL  STOP' 
GOTO  225 
END  IF 

YLDl  -  halfPI'SIGMAX/ 

DSQRT (2. /DCOS(halfPx*SIGMAX/TYLD) -2.) 

if  ( (a-aol)  .It.  wol)  then 

CODX  -  CODMAX' (1 (A-AOL) /WOL) **2 
IF  (CODX. LT. 0.0)  CODX  -  0.0 
CHECK  -  DSQRT (2. 'E*YLD1*C0DX/BETAI ) 
end  i  f 

IF  ((  (A-AOL)  .  ge  .  WOL)  .  or .  (Iimax  .gt .  eheclt)  )  then 
KMAXOL  -  KMAX 
KMINOL  -  KMIN 


o  n  n  o  o  o 


A.IO 


WOL  -  PI/8. ‘BETA* (KMAX0L/YLD1)**2 
CODMAX  -  BETAl*KMAXOL**2/2./E/yi,Dl 
AOL  -  A 
end  if 
end  if 

IF  (KMIN.le.KHINOL)  then 
KMINOL  -  KMIN 
SR  -  KHINOL/KMAXOL 

XR  -  0. 46733+0. 29'I01*SR+0.23866*SR**2 
XALPHA  -  l.-(l.-DSORT(BETA))*<I.-SR) 

XFR  -  l.-(0.7S6-0.912*FR+0.156*FR**2)* (l.-SR) 

EFFKC  -  XR*XALPHA*XFR*KMAXOL  -  0.225*KMINOL 
IF  (SR.LT.O.)  EFFKC  -  EFFKR*KMAXOL*DEXP (0. OB'SR) 
end  if 

REFK  -  AMP*SMAX1*C0RFAC 
IF  (RKFK  .ge.  KC)  GOTO  210 
C 

C  COMPUTE  EFFECTIVE  KMIN  AND  DELTA  K 

C 

IF  (KMAX  .le.  0.)  THEN 
c  . . .da/dn  -  0 

X2  -  A 
ELSE 

KMXEFF  ■-  EFFKC+0.225*KMAX 
IF  (KMXEFF.lt. KMAX)  KMXEFF  -  KMAX 
IF  (KMXEFF. LT. 0. )  KMXEFF  -  0. 

IF  (KMAX.GE.KC)  THEN 

PRINT  '.'KMAX  EXCEEDED  KC'.KMAX 
GOTO  210 
END  IF 

KMNEFF  -  EFFKC+0.225*KHIN 
IF  (KMNEFF. LT.O.)  KMNEFF  -  0. 

IF  (KMNEFF.lt. KMIN)  KMNEFF  -  KMIN 
DKEFF  -  KMXEFF-KMNEFF 

IF  KMAX  IS  LESS  THAN  KMNEFF  GROWTH  IS  ZERO. . . 

IF  (DKEFF.gt.O.)  THEN 

COMPUTE  DA/DN  AND  CRACK  GROWTH. . . 

DKEFF  -  DLOGIO (DKEFF) 

CALL  ILU (DELKB, DADNB, NDADN, DKEFF, DADN) 
cc  old:  DADN  -  10. ••DADN'KC/ (KC-KMAX) 

dadn  -  10.**dadn 
X2  "A 

A  =  A  +  DADN'CYCLES 

IF  (A  .ge.  W/2)  THEN 

REFK  -  AMP*SMAX1*C0RFAC 

PRINT  *,'THE  CRACK  LENGTH  IS  EXCEEDING  THE  WIDTH' 
PRINT  *, 'HOURS  -  HOURS 
PRINT  *. 'LI/2  -  '.A 

PRINT  '.'REFK  -  ' , REFK 
GOTO  250 
END  IF 

END  IF 
END  IF 

200  CONTINUE 


C 

C  CHECK  THAT  KC  HAS  NOT  BEEN  EXCEEDED... 

IF  (REFK.GT. . 95*KC)  THEN 
CALL  KCOR(a,w,CORFAC) 

REFK  -  AMP*SMAX1’>C0RFAC*DSQRT(PI*A) 
END  IF 

IF  (REFK  .It.  KC)  GOTO  225 

210  XONE  -  X2 
DO  J  -  3,7 

o  Initialize... 

XINCR  -  l./(10.»*J) 

REFKIA  -  Kc  -  1. 

II  -  I 


non 


A.ll 


DO  while  ( (REFKIA.le.Kc)  .and.  (II.lc.lOO)) 

A  -  XONE+II*XINCR 

if  (a  .ge.  w/2)  then 

REFK  -  AMP*SMAXl*CORFAC 

PRINT  *,'THE  CRACK  LENGTH  IS  EXCEEDING  THE  WIDTH' 
PRINT  *, 'HOURS  -  HOURS 
PRINT  *, 'LI/2  -  '.A 

PRINT  *, 'REFK  -  '.REFK 
GOTO  250 
end  if 

CALL  KCOR(a,w,CORFAC) 

CORFAC  «  CORFAC*DSORT(PI*A) 

REFKIA  -  AMP'SMAXl 'CORFAC 
II  -  II+l 
END  DO 

IF  <II.GE.100)  THEN 

PRINT  *, 'PROBLEMS  WITH  ACRR  ITERATIONS' 

STOP 

r  T  CP 

XONE  -  A-XINCR 
END  IF 
END  DO 

IF  (M.EQ.l  .AND.  N.EQ.l)  THEN 
PRINT  *, 'FAILED  IMMEDIATELY' 

IFALl  -  1 
GOTO  250 
END  IF 

HOURS  »  BHRS* (M-1+ ( (YCYC-CYCLES)/XCYC) ) 

REFK  =  REFKIA 

WRITE  (6,320) 

WRITE  (6,350)  HOURS, A, REFK 
GOTO  250 
(OUTPUT  RESULTS) 


225  HOURS  •  BHRS'M 

IF  (HOURS. ge.PHOURS)  then 
PHOURS  -  PHOURS+PRTHRS 
WRITE  (6,350)  HOURS, A, REFK 
KK  -  KK+1 

HI(KK)  =■  HOURS 
AI (KK)  =  A 
end  i  f 

REWIND  3 
YCYC  =  0. 

END  DO 

250  IF  (IFALl. ne.l)  then 
KK  =  KK+1 

HI  (KK)  =  HOURS 
AI (KK)  =  A 
end  i  f 


CLOSE  (3, STATUS- 'DELETE ' ) 


888 


write  (• , ' ( /lx, a, S) ' )  'Enter  the  name  of  your  PREDICTION  file:  ' 
read  ( * ,  '  (  a  ) ' )  ipred 

open  (unit=l,  file-ipred,  status- ■ new' ,  err-888) 


write 

(1, 

ax,a/)  ') 

*  INPUT  DATA’ 

write 

(1, 

<lx,a/) ') 

titles 

write 

(1, 

{lx, a/lx . a/ ) 

' ) 

’No  rORMAN-transformation  of  dadn-corve 
’and  YLDl  calculated  outside  if-block' 

write 

(1, 

{lx, a,  lx, a 

)  ' ) 

'Analysis  data:  file  iafid 

write 

(1, 

{lx, a,  F7.3 

)  ') 

•  ALPHA  alpha 

write 

(1, 

(lx,a,E13.5 

)  ’) 

*  Ain  «* ,  din 

write 

(1, 

(lx,a,E13.5 

)  ') 

*  W  *• ,  w 

write 

(1, 

(ix, a, El  3 . 5 

)  ' ) 

*  DLS  «»,  DLS 

write 

<1, 

(lx, a,  F7.3/)') 

•  AMP  anp 

write 

(1, 

(lx, a,  lx, a  )  ' ) 

•Material 

data : 

file 

-',  mat fll 

write 

(1, 

(lx, a, Ell. 3  ) ' ) 

1 

E 

E 

write 

(1, 

(lx, a,  F7.3  ) ' ) 

1 

Xnu 

-•,  xnu 

write 

(1, 

(lx, a, Ell. 3  ) ' ) 

t 

Tyld 

-',  tyld 

write 

(1. 

(lx, a.  F7.3  )  ') 

1 

FR 

-•,  FR 

write 

(1. 

(lx, a, Ell. 3/) •) 

« 

Kc 

-',  KC 

write 

(1, 

(lx, a,  lx, a  )  ■) 

•Spectrum 

data: 

file 

-',  isfld 

write 

(1, 

(lx, a,  17  )•) 

« 

Nlevs 

-',  nlevs 

write 

(1, 

(lx, a, Ell. 3  ) ') 

* 

Bhrs 

-',  bhrs 

write 

(1, 

(lx, a, Ell. 3/) ■ ) 

1 

Smaxl 

smaxl 

write 

(1, 

(/lx, a/) ') 

1 

OUTPUT  DATA 

IF  (IFALl.ne.l)  then 

WRITE  (1.997)  H.ITCYC.N, HOURS 
WRITE  (1,996)  A.REFK 
WRITE  (1,995)  SIGMAX.SIGMIN, CYCLES 
end  If 

write  (1 , '  (/lx, a) ' )  '  Hours  A' 

do  i  -  1,KK 

write  (1,  *  (lx,F19.3,F15.6)  •)  HI(i),AI(i) 
end  do 

CLOSE  (1) 

STOP  1776 

280  FORMAT (  /  ’  CONTACT  STRESS  MODEL-III  ANALYSIS'  /) 

320  FORMAT (/2X, 'STRESS  INTENSITY  EXCEEDED  IN  A  DIRECTION  :'/) 

339  FORMAT(/'  BEGIN  SPECTRUM  CRACK  GROWTH  ANALYSIS,  A  -',E12.5 
+//T8, 'HOURS', T22, 'A',T29, '  KREFERENCE' /) 

350  FORMATC  ' ,  Fll .  0,  FIS  .  9,  F12 . 1 ) 

997  FORMAT(  5X, 'BLOCKS  - ' , 18, /5X, ' TOT . CYC  -', 18, /5X, ' LAYER  -'18, 
♦  /5X, 'HOURS  -',F9.0) 

996  FORMAT (SX. 'A  - ' , F9 . 6, /5X, ' KREFA  - ' , ell . 4, /5X) 

995  FORMAT (SX, 'SIGMAX  - ' , ell . 4 , /5X, ' SIGMIN  -', ell . 4, /5X, 'CYCLES 
+  ,F10.0,/) 


END 


SUBROUTINE  KCOR ( a , w , CORFAC ) 

IMPLICIT  REAL*8  (A-H,0-Z) 

DATA  PI/3.1415926536/ 

CORFAC  -  DSQBT(l./dcos(pi*a/w) ) 

RETURN 

END 

SUBROUTINE  TLU (X, Y, N, XVAL, YVAL) 

IMPLICIT  REAL’S  (A-H,0-Z) 

DIMENSION  X (N) , Y (N) 

IF(X (N) .GE.XVAL)  GO  TO  10 
I  -  N 
GO  TO  30 

10  DO  20  I  =  1,N 

IF(X(I)  .GE.XVAL)  GO  TO  30 
20  CONTINUE 
30  IF  (I.EO.l)  1-2 

YVAL  -  (Y  (I) -Y(I-l)  )  /  (X(I) -X  (I-l)  )  •  (XVAL-X(I-l)  )+Y(I-l) 

RETURN 

END 


Examole  problem 
A 

.1 


0 

0 

dadn.dac 

•too 

100. 


Table  A.l  Example  analysis  file,  named  example. anl. 


Material's  curve  for  At-UISG-TSSl  in  [N,ml 

6.90el0  .3 

457. e6  411. e6 

9.9Qe7  14 

1454550. 

l.OOOOOOE-011 

1711900. 

4.000000E-011 

2201980. 

1.993473E-010 

3097420. 

9.201416E-010 

2.69774E+007 

2.86472SE-006 

3.30770E+007 

8.310194E-006 

3.70996E+007 

2.461269E-005 

3.99125E+007 

7.3B9003E-005 

4.19771E+007 

2.238263E-004 

4.42058E+007 

1.193590E-003 

4.57e76E+007 

6.425438E-003 

4.69660E+007 

3.48196SE-002 

4.78767E+00'7 

1.896133E-001 

4.86000E+007 

1.0 

Table  A.2  Example  material  file,  named  dadn.dat. 


Example 

spectrum  file 

(1) 

5 

100. 

1.00 

-0.20 

1.00 

0.05 

1000. 

0.40 

-0.10 

5000. 

0.60 

0.00 

3000. 

1.00 

0.40 

2000. 

0.50 

-0.20 

3500. 

0 

Table  A. 3a  First  example  spectrum  file,  named  exampll.seq. 


Example 

spectrum 

file  (2) 

5 

100, 

2.00 

-0.20 

1.00 

0.C5 

1000. 

0.40 

-0.10 

5000. 

0.60 

0.00 

3000. 

1.00 

0.40 

2000. 

0.50 

-0.20 

3500, 

0 

Tabic  A.3b 


Second  example  spectrum  file,  named  exampll.seq. 


Modified  version  of  the  McAir  crack  growth  program  CNTKM8. 

No  Forman-transf ormation  of  the  da/dn-curve 

Marcel  Bos,  5/3/91. 

WARNING:  In  case  of  pre-cycled  specimens  it  does  not  matter 

whether  Sl-units  or  Imperial  units  are  used,  as  long 
the  input  is  consistent.  In  case  of  saw-cut  specimens 
(line  3  of  spectrum  file:  0,0)  only  Sl-units  should 
be  used  however,  although  the  error  that  is  introduced 
when  using  imperial  units  most  likely  is  very  small. 


Input  ANALYSIS  file:  example.anl 
Input  SPECTRUM  file:  exampll .seq 

CONTACT  STRESS  MODEL-lII  ANALYSIS 
Example  problem 

Input  ALPHA:  0 
AINITIAL:  .02 
DLS:  60e6 
AMP  FACTOR:  1 

BEGIN  SPECTRUM  CRACK  GROWTH  ANALYSIS,  A  -  .20000E-01 


HOURS 

A 

KREFERENCE 

100. 

.021410967 

17438515.0 

200. 

.023101030 

18484567.7 

300. 

.025263967 

19895103.3 

400. 

.028267137 

22023761.4 

500. 

.033275489 

26150467.9 

STRESS  INTENSITY  EXCEEDED  IN  A  DIRECTION  : 
576.  .043440339  49000029.5 

Enter  the  name  of  your  PREDICTION  file:  predl 
Return  code  1776 


Tabic  A.4a  Interactive  input  of  first  example.  User's  input  written  in  italic. 


Modified  version  of  the  HcAlr  crack  growth  program  CNTKM8. 

No  Forman-transformation  of  the  da/dn-curve 

Marcel  Bos,  5/3/91. 

WARNING;  In  case  of  pre-cycled  specimens  it  does  not  matter 

whether  Sl-units  or  Imperial  units  are  used,  as  long 
the  input  is  consistent.  In  case  of  saw-cut  specimens 
(line  3  of  spectrum  file:  0,0)  only  Sl-units  should 
be  used  however,  although  the  error  that  is  introduced 
when  using  imperial  units  most  likely  is  very  small. 


Input  ANALYSIS  file:  example.onl 
Input  SPECTRUM  file:  exampG.scq 

CONTACT  STRESS  MODEL-IIl  ANALYSIS 
Example  problem 

Input  ALPHA;  0 
AINITIAL:  .02 
DLS:  60e6 
AMP  FACTOR:  _5 

BEGIN  SPECTRUM  CRACK  GROWTH  ANALYSIS,  A  -  .20000E-01 

HOURS  A  KREFERENCE 

100.  .020103570  16779813.0 

200.  .020212387  16846233.6 

300.  .020325504  16915458.4 

400.  .020443242  16987709.7 

500.  .020565959  17063236.6 

600.  .022061815  17836254.6 

700.  .023919250  19007745.9 

800.  .026362594  20648494.6 

900.  .029945952  23318507.4 

1000.  .037060907  29836113.9 

STRESS  INTENSITY  EXCEEDED  IN  A  DIRECTION  : 

1007.  .043440407  49000320.5 

Enter  the  name  of  your  PREDICTION  file:  pred2 
Return  code  1'776 


TftWe  A.4b  Interactive  input  of  second  example.  User’s  input  written  in  italic. 


INPUT 

DATA 

Example  problem 

No  rORMAN-transformatlon  of  dadn-curve 

and  YLDl  calculated  outside  If-block 

Analysis  data: 

file  -  example. anl 

ALPHA  -  .000 

Ain  -  .20000E-01 

W  -  .lOOOOE+00 

DLS  -  .60000E+08 

AMP  -  1 . 000 

Material  data: 

file  -  dadn.dat 

E  -  .690E+11 

Xnu  -  .300 

Tyld  -  .457E+09 

FR  -  .899 

Kc  -  .490E+08 

Spectrum  data: 

file  «  exampll.seq 

Nlevs  ”  5 

Bhrs  -  .lOOE+Oa 

Smaxl  -  .600E+08 

OUTPUT 

DATA 

BLOCKS  - 

6 

TOT. CYC  • 

87000 

LAYER 

5  *) 

HOURS 

576. 

A  • 

043440 

KREFA 

4900E+08 

SIGMAX  - 

3000E+08 

SIGMIN  -  - 

1200E+08 

CYCLES  - 

3500. 

Hour  j 

A 

.000 

.020000 

100.000 

.021411 

200.000 

.023101 

300.000 

.025264 

400.000 

.028267 

500.000 

.033275 

575.862 

.043440 

'Layer'  is  (he  McAir  dcnoution  for 'block'  as  used  in  this  report,  while 


'block'  in  McAir  jargon  deirates  'sequence'  or 'program'. 


Table  A.5a  Prediction  file  of  first  example,  named  predl. 


INPUT 

DATA 

Example  problem 

No  FORKAN-tranaformation  of  dadn-curve 

end  YLDl  calculated  outside  if-blocX 

Analysis  data: 

file  -  example. anl 

ALPHA  -  .000 

Ain  -  .20000E-01 

W  -  .lOOOOE+00 

DLS  -  .60000E+08 

AMP  -  .500 

Material  data: 

file  -  dadn.dat 

E  -  .690E+11 

Xnu  -  .300 

Tyld  -  .457E+09 

FR  -  .899 

Kc  -  .490E+08 

Spectrum  data: 

file  “  exampl2.seq 

Nlevs  -  5 

Bhrs  -  .lOOE+03 

Smaxl  -  .120E+09 

OUTPUT 

DATA 

BLOCKS  * 

11 

TOT. CYC  -  151000 

L/.YER 

2 

HOURS 

1007. 

A 

043440 

KREFA  - 

4900E+08 

SIGHAX  « 

2400E+08 

SIGMIN  =■ 

6000E+07 

CYCLES  - 

5000. 

Hours 

A 

,000 

.020000 

100.000 

.020104 

200,000 

.020212 

300,000 

.020326 

400.000 

.020443 

500,000 

.020566 

600,000 

.022062 

700.000 

.023919 

800,000 

.026363 

900.000 

.029946 

1000.000 

.037061 

1006.897 

.043440 

Table  A.5b  Prediction  file  of  second  example,  named  predl. 
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Flow  diagram  of  computer  program  CNT8SHR. 
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Figure  A.l 


(continued) 
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APPENDIX  B  COMPUTER  PROGRAM  ARLECC 

In  this  appendix  the  source  code  of  computer  program  ARLECC  is  listed.  This  program 
incorporates  the  crack  closure  model  as  proposed  in  chapter  4  of  this  report.  It  is  used  for  the  calculation 
of  the  crack  growth  life  of  CCT  specimens  with  through  cracks.  It  has  been  written  in  Fortran  77,  and  it 
can  be  run  on  a  personal  computer  with  a  mathematical  co-processor.  Although  the  same  spectrum  files 
can  be  used  as  for  program  CNT8SHR,  even  small  jobs  with  simple  blocked  spectra  may  take  more  than 
an  hour,  depending  on  the  applied  stress  levels.  This  is  due  to  the  fact  that  ARLECC  treats  cycles 
individually,  allowing  the  crack  to  grow  each  cycle,  whereas  CNT8SHR  operates  on  the  load  blocks  as  a 
whole,  although  it  is  possible  to  perform  a  cycle-by-cycle  calculation  with  the  latter  program  by 
assigning  only  one  load  cycle  per  block.  Unlike  CNT8SHR,  computer  program  ARLECC  always 
assumes  that  the  initial  crack  has  experienced  prior  growth,  so  that  residual  plastic  zones  are  present 
around  the  crack  tip.  Computer  program  ARLECC  can  also  be  used  to  perform  'no-load-interaction 
analyses',  in  which  case  only  the  stress  ratio  effect  is  accounted  for,  and  previous  load  history  is  ignored. 
Any  consistent  set  of  units  can  be  used  for  the  input  data. 

B.l  Program  input 

The  input  consists  of  two  input  files  and  an  interactive  part.  The  two  input  files  arc  the  material 
file  and  the  spectrum  file.  They  arc  identical  to  the  ones  needed  for  computer  program  CNT8SHR, 
Details  can  be  found  in  appendix  A.  Again,  a  crack  growth  curve  for  a  suess  ratio  of  zero  is  expected. 
The  interactive  part  is  very  straightforward.  Successively  the  program  asks  for:  The  name  of  the 
spectrum  file,  the  name  of  the  material  file,  the  value  of  retardation  parameter  a,  the  initial  half  crack 
length,  the  specimen  width,  the  design  limit  stress,  the  amplification  factor,  and  the  maximum  number 
of  programs  to  be  analyzed.  All  these  parameters,  including  a,  are  similar  to  the  ones  described  in 
appendix  A.  To  perform  a  'no-load-interaction'  analysis  a  value  of  -999  must  be  input  for  a.  In  all  other 
cases  a  value  of  between  -1  and  +1  must  be  used. 

Table  B.l  shows  what  the  screen  looks  like  when  running  the  first  e.xample  of  appendix  A.  Again 
the  data  files  of  tables  A.2  and  A.3a  arc  used.  The  output  file  pred]  is  given  in  table  B.2. 

B.2  Listing  of  ARLECC 

Most  of  the  variable  names  arc  identical  to  the  ones  used  in  computer  program  CNT8SHR.  The 
source  code  is  amply  provided  with  commentary  lines.  Figure  B.l  shows  the  How  diagram  of  the 
program. 
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progrcm  ARLCCC 

ccccccccccccccccccccccccceccccccccccccccccccccccceccccccccccccccccccccc 


c  c 
c  Program  to  calculate  crack  growth  under  VA-loading  in  c 
c  CCT-specimens.  Based  on  equations  duo  to  Schijve  and  McAir.  c 
c  c 
c  Marcel  Bos,  17/4/91  c 
c  c 


cccccccccccccccccccccccocccccccoccccccccccccccccccccccccccccccccccccccc 

implicit  none 

integer  Icount, 11, j j, kk, nblks, ndadn, nlevs 
real  *8  a,ai(500),ain,  alpha,  amp,  aol,beta,  bhrs, 

+  check, corfac, cycles, eyld, dummy, da, 

+  dadn(50)  ,delk  (50) ,  dk, dkef f , dls,  fO,  fr,hl  (500) , 

+  hours, kc, kmax, kmaxol,kmin,kmlnol, kmnefl, 

+  pi, ref k, sigmax, slgmin, smax, smin, 

+  ty Id, totcyc, xcyc, w, wol , yld, yldol 

logical  failed, unrctarded 

character  ipred"20, i sf ld*20, mat f 11*20, title*72 
data  pi  /3. 1415926536/ 

call  HEADER 

write  (*, • (/12x,a/12x,a/)  •) 

+  'ARL  craclt  growth  program  for  CCT  specimens  under  VA  loading.', 

+  'Make  sure  that  the  da/dn-dK  data  are  for  stress  ratio  R  -  0.' 
write  (*, '  (24x,a/) ' )  '  Marcel  Bos,  17/4/91' 

c 

c  read  data . . . 
c 

22  write  (*, '  (lx, a, $) ' )  'Input  SPECTRUM  file:  ' 
read  (* ,  '  (  a  )  ' )  isf  id 

open  (unit-4,  flle-isfid,  err-22,  status-'old' ) 

write  (*, •  (lx, a, $) ' )  'Input  MATERIAL  file:  • 
read  (*, •  (a) ' )  matfll 

33  write  (*, ' (/lx, a, $) ' )  ‘Input  ALPHA:  ' 
read  (*,  *  )  alpha 

if  (alpha  .eq.  -999.)  then 
unretarded  «  .true, 
else 

if  (  (alpha. ge. -1 . 0)  .and.  (alpha. le. 1 .0) )  then 
unretarded  -  .false, 
else 

go  to  33 


end 

if 

end  1  f 

write 

*,  •  (  lx,a, S)  •)  '  AINITIAL:  ' 

read 

*,  *  ) 

ain 

write 

*, ' (  lx,a, S)  ' )  • 

W:  • 

read  i 

*,  *  ) 

w 

write  1 

*,  '  (  lx,a,S)  ')  • 

DLS:  • 

read  i 

*,  •  ) 

dls 

write  1 

*, ' (  lx,a, $) ’ )  •  AMP 

FACTOR:  ’ 

read  i 

*,  *  I 

amp 

write  1 

*, '  (  lx,a,$)  ')  ' 

NBLKS:  ' 

nblks  - 

500 

write  1 

*,  '  (i3/)  ' )  nblks 

if  (dls 

.le.  0.)  stop  9991 

open  (unit-1,  file-matfil. 

status»*old 

read  (1, '  (a)  ' )  title 

read  (1,  *  )  dummy, dummy , tyid, cyld 

read  (1,  *  )  Kc,ndadn 


if  (Kc  .le.  0.0)  stop  9992 

fr  ”  cyld/tyld 
fr  -  fr/1.12 

...this  is  done  to  facilitate  Schijve's  equation  in 
subroutine  effkmln,  which  is  valid  for  2024  and  7075, 
these  materials  having  a  fr  of  approx.  1.12 

to  -  1.  -  0.45*(0.244  +  0.912*fr  -  0.156*fr*fr) 
c  ...f0*dK  yields  the  effective  stress  intensity  for  R  -  0, 

c  see  also  subroutine  etfkmin 

do  11  -  l.ndadn 
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r««d  (1, 
dadndl) 
delk(li) 
end  do 


delk  (11) ,dedn  (11) 
dloglO(  dadn(ll)  ) 

dlo(]10(  fO«del)c(il)  ) 


close  (1) 


e 

c 

c 


reed  stress  exceedance  date... 

write  (*,'(lx,e)')  'Initializing  temporary  stress' 

write  (*,'(lx,a)')  'spectrum  file  (this  may  talce  a  while)...' 

open  (unlt-3,  flle-'temp',  atatus-'new' ,  form-'unformatted' ) 

road  (a,  '  (a) ' )  title 
read  (4,  *  )  nlevs.bhrs 

read  (4,  *  )  smax.smln 

smax  •  smax*dls 
smin  •  smln*dls 
xcyc  -  0. 

do  11  •  l.nlcvs 

read  (4,*)  slgmax, slgmln. cycles 
xcyc  -  xcyc+cycles 
slgmax  •  slgmax*dls 
sigmin  -  sigmin*dls 
if  (slgmax. gt . smax)  smax  -  slgmax 
if  (sigmin. It . smin)  smln  ■  sigmin 
write  (3)  slgmax, sigmin, cycles 
end  do 

If  (smax  .le.  0.)  stop  9994 

close  (4) 
close  (3) 


if  (.not. unretarded)  then 

beta  -  0.333333  +  0.5*alpha  <•  0.166667*alpha*alpha 
...alpha  ■  -1  is  minimum  value  for  alpha;  this  approximates 
unretarded  case  but  is  not  the  same  however; 
input  alpha  -  -999  for  that  case! 

alpha  «  0  - >  plane  strain,  beta  »  1/3 

alpha  ■  1  - >  plane  stress,  beta  »  1 

el  se 

beta  •  0. 
end  if 

initialize. .  . 


totcyc 

icount 

a 

al  (1) 
hi  (1) 
failed 


0. 

1 

aln 

a 

0. 

. false. 


corfac  -  dsqrt (1. /dcos (pi 'a/w) ) 
corfac  -  corfac'dsqrt (pl*a) 
refk  -  amp'corfac'smax 
if  (ref!«  .ge.  )cc)  failed  -  .true, 
if  (failed)  go  to  999 


assuming  that  the  specimen  already 
experienced  crac)c  growth... 
kmaxol  -  corfac*3max 
kmlnol  -  cortac'smin 


yldol  -  . 5*pi*smax/dsqrt (2. /dcos (0.5*pi*smax/tyld) -2.) 
aol  -  a 

wol  -  (beta/pl )  *  (limaxol/yldol)  ••2. 
calculate  min.  eff.  stress  intensity  )cmneff... 


call  EFFKMIN  ()cmlnol,kmaxol,  fr,  Icmnef  f ) 


start  crack  growth... 

write  ( ' ,  '  (/lx.  a/ )  ' )  'Et.irt  crack  growth  calculation...' 

write  ( • ,  '  ( 6x, as, 1  Ox, al , 1  Ox, a4 ) ■ )  ' hours ' ,  ' a ' ,  ' ref k ' 

open  (unlt-3,  flie-'temp',  status- ' old’ ,  form-'unformatted') 
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do  11  •  l,nbl)ca 
do  jj  •  l.nlevs 

read  (unlt-3)  slgmax, slgmln, cycles 
if  (sigmax  .le.  slgmln)  stop  9998 
If  (sigmax  .ge.  tyld)  stop  9999 

do  kk  ■  1. nint (cycles) 

totcyc  -  totcyc  +  1. 

kmax  -  cor fac* sigmax 

kmin  "  corfac'slgmin 

dk  -  kmax  -  kmin 

if  (unretarded)  then 

call  ErFKMIN (kmin, kmax, fr,kmneff) 

...kmneff  of  the  current  cycle  is  used 
else 

...kmneff  of  the  previous  cycle  is  used 
end  if 

dkcff  -  kmax  -  kmneff 

if  ((dkeff  .gt.  dk)  .and.  (kk  .ne.  1))  dkeff  -  dk 

if  (dkeff  .gt.  0.001'dK)  then 

dkeff  -  dloglO (dkeff) 

call  TLU (delk,dadn,ndadn, dkeff, da) 

da  -  10. "‘da 
a  -  a  *  da 

if  (a  .ge.  w/2)  then 
failed  •  .true, 
else 

corfac  -  dsqrt (1 . /dcos (pi*a/w) ) 
corfac  -  corfac'dsqrt  (pi*a) 
refk  -  amp*corfac*smax 
if  (refk  .ge.  kc)  failed  •  .t  i"ie. 
end  i  f 

if  (failed)  go  to  999 
end  if 

now  check  whether  current  load  sigmax  is  overload... 

if  (. not .unretarded)  then 

if  (sigmax  .ge.  0.001)  then 

yld  =  0. 5*pi*sigmax/ 

dsqrt (2 . /dcos (0 . 5*pi * sigmax/ty Id)  -  2.) 

if  (a  .It.  aol+wol)  then 

check  -  kmaxol *dsqrt (1 . - (a-aol) /wol ) 'yld/yldol 
else 

check  -  0. 
end  i  f 

if  (kmax  .gt.  check)  then 
kmaxol  •  kmax 
kmlnol  ■  kmin 
yldol  =■  yld 
aol  'a 

wol  -  (beta/pi) * (kmaxol/yldol) **2. 

end  if 
end  1  f 

...and  check  whether  sigmin  is  underload... 

(remember  that  kminol  has  already  been  set  to  kmin  in 
case  of  sigmax  being  an  overload) 

if  (kmin  .le.  kminol  +  .OOl'dK)  then 

kminol  -  kmin 

calculate  min.  eff.  stress  Intensity  kmneff... 

call  EFFKMIN (kminol , kmaxol,  f r,  kmneff) 
end  1  f 


end  I 
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end  do 
end  do 

hours  «  bhrs*H 

icount  -  ieount  +  1 
•l(lcount)  -  a 
hi<lcount)  -  hours 

write  (*, ’ (lx, flO.2, 2x, f9. 6, 2x,el2.4) • )  hours, a, refX 
rewind  (3) 
end  do 

999  if  (failed)  then 
c  ...specimen  has  failed 

write  (*, ' (/lx, a/) ' )  '...specimen  has  failed' 
icount  •  Icount  +  1 
ai  (icount;  -  a 

hi  (icount)  -  bhrs*  (totcyc/xcyc) 
else 

c  ...nbl)cs  was  exceeded 

write  (* ,  '  (/lx,  a/)  ' )  '...maximum  number  of  bloc)(5  exceeded' 
end  if 

close  (3, statusa'delete' ) 


c  output . . , 
c 

94  write  (• , ' (lx, a, S) ' )  'Enter  the  name  of  your  PREDICTION  file:  • 
read  (* ,  '  (  a  ) • )  ipred 

open  (unit«l,  file-ipred,  status-'new' ,  err-94) 

write  (1, '  (lx,a/) ')  'ARL  Crach  Growth  Model ... ' 

write  (1, '  (lx,a/) ■)  •  INPOT  DATA' 

if  (unretarded)  then 
write  (1, • (lx, a) ' ) 

+  '  ALPHA  •  N.A.  (unretarded)' 

else 

write  (1, ' (lx, a,  F7.3  ) ') 

+  '  ALPHA  alpha 

end  i  t 


write 

(1, 

•  (lx,a,E13.5  )  ') 

• 

Ain 

ain 

write 

(1, 

'  (lx, a, E13 . 5  )  ' ) 

« 

W  , 

w 

write 

(1, 

(lx,a,E13.5  )  ') 

1 

DLS 

DLS 

write 

(1, 

(lx, a,  F7.3/) ') 

1 

AMP  »*, 

amp 

write 

(1, 

(lx, a,  lx, a  )  ' ) 

'Material  data: 

file 

matfil 

write 

(1, 

(lx, a, Ell. 3  )■) 

• 

Tyld 

tyld 

write 

(1, 

(lx, a,  FT. 3  )  •) 

• 

FR  »•, 

FR 

write 

11. 

(lx, a, Ell. 3/)  ') 

• 

Kc 

KC 

write 

(1, 

(lx, a,  lx, a  )  ' ) 

’Spectrum  data: 

file  -* , 

isf  id 

write 

<1, 

(lx, a,  17  )') 

Nlevs 

nlevs 

write 

(1, 

(lx, a, Ell. 3  )  •) 

s 

Bhrs  ■’ , 

bhrs 

write 

(1, 

(lx, a, Ell. 3/) •) 

• 

Smax  - ’ , 

smax 

write 

(1, 

(/lx, a/)  •) 

'  OUTPUT 

DATA’ 

write 

(1, 

(/lx, a)  ' )  ' 

Hours 

A' 

do  ii 

-  1 

Icount 

write  (1,  ■  (lx,  fl<1.3,  flS.6)  • )  hiiii) ,  ai  (ii) 
end  do 


close  (1) 

stop 

END 


subroutine  TLO (x, y , n, xval , yval ) 

Implicit  real*8  (a-h,o-r) 
dimension  x(50),y(50) 

i f (x (n) . ge . xval )  go  to  IQ 
i  •  n 
go  to  30 

10  do  20  1  -  l,n 
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if (X  (i) .ge-xvaX)  go  to  30 
20  continue 
30  if  (i.eq.l)  i  -  2 

yval  -  (y(i)-y(i-i))/(x(i)-x(i-l))*(xval-x{i-l))+y(i-l) 

return 

END 


subroutine  EFFKHIN (kminol, Wmaxol, fr, kmnef f ) 
implicit  none 

real *8  kminol, kmaxol, fr, kmnef f ,  SR, U,XR,XFR 

SR  -  Xminol/Xmaxol 
U  -  0.55  +  0.33*SR  +  0.12*SR*SR 
XR  -  1.0  -  U« (1.-  SR) 

XFR  -  1.0  -  (0.756  -  0.912*FR  +  0 . 156*FR*FR) • <1 . -  SR) 

c  . U  -  dkeff/d)c  is  an  equation  due  to  Sc)iljve 

c  while  the  equation  for  XFR  has  been  obtained  from  McAlr. 
c  To  tune  the  combination  of  these  equations  FR  is  divided  by 
c  a  factor  1.12  in  the  main  program  since  Schijve's  equation  is 

c  valid  for  2024  and  7075,  having  a  FR  of  1.12  approx. 

)tmneff  ■  XR*XFR*)cmaxol 

if  Ocmneff  .It.  0.)  )cmneff  ”  0. 

if  Ocmneff  .It.  )cminol)  Icmneff  ■■  jcminol 

return 

END 
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ARL  crack 

Enhanced  Crack  Closure  model  ARLECC 

growth  program  for  CCT  specimens  under 

VA  loading 

Hake  sure 

that  the  da/dn-dK  data  are  for  stress 

ratio  R  •  0 

Marcel  Bos,  17/4/91 


Input  SPECTRUM  file:  eXOmpll .SCq 
Input  MATERIAL  file;  dodn.dot 

Input  ALPHA:  0 
AINITIAL:  .02 
W;  ./ 

OLS:  60e6 
AMP  FACTOR:  1 
N3LKS:  500 

Inltializinq  temporary  stress 
spectrum  file  (this  may  take  a  while)... 

Start  crack  growth  calculation... 


hours 

a 

refk 

100.00 

.020908 

.1728E+08 

200.00 

.021942 

.1793E*08 

300.00 

.023141 

.1872E+08 

400.00 

.024566 

.1969E+08 

500.00 

.026326 

.2097E*09 

COO. 00 

.028635 

.2282E+08 

700.00 

.032063 

.2606E+08 

800.00 

.041695 

.4276E+0B 

...specimen  has  failed 


Enter  the  name  of  your  PREDICTION  file:  predi 
Stop  -  yroqram  terfnln^icd. 


Tabic  C.l  Inicniciivc  input  of  first  example.  User's  input  written  in  italic. 


ARL  Crack  Growth  Model 


INPUT  DATA 


Material  data: 


Spectrum  data: 


ALPHA 

.000 

Ain 

.20000E-01 

w 

.lOOOOE+00 

DLS 

.SOCOOE-t-OS 

AMP 

1.000 

file 

dadn.dat 

Tyld 

.457E+09 

PR 

.803 

Kc 

.490E+08 

file 

exampll.seq 

Nlevs 

5 

Bhre 

.lOOE+03 

Smax 

.eooE+oe 

OUTPUT  DATA 


Hours 
.000 
100.000 
200.000 
300.000 
400.000 
SOO.OOO 
600.000 
700.000 
BOO. 000 
800.007 


A 

.020000 

.020908 

.021942 

.023141 

.024566 

.026326 

.028635 

.032063 

.041695 

.056124 


Table  B.2  Prediction  file  of  first  example,  named  predl. 


Figure  B,  I  Flow  dinf;ram  of  computer  program  ARLECC. 


fjKurcB.I  (continued) 


DISTRIBUTION 


AUSTRALIA 

l>gpannignt-pfPgfgng.s 

Defence  Central 

Chief  Defence  Scientist  ) 

AS,  Science  Corporate  Management  )  shared  copy 
FAS  Science  Policy  ) 

Director,  Departmental  Publications 

Counsellor,  Defence  Science,  London  (Doc  Data  sheet  only) 

Counsellor,  Defence  Science,  Washington  (Doc  Data  sheet  only) 
Scientific  Adviser,  Defence  Central 
OIC  TRS,  Defence  Central  Library 
Document  Exchange  Centre,  DSTIC  (8  copies) 

Defence  Intelligence  Organisation 

Librarian  H  Block,  Victoria  Barracks,  Melbourne  (Doc  Data  sheet  only) 

Aeronautical  Research  Laboratory 
Director 
Library 

Chief  -  Aircraft  Structures  Division 
Divisional  File  -  Aircraft  Structures 
Author:  M.J.  Bos 

C. K.  Rider 
J.G.  Sparrow 

J.M.  Finney  (2  copies) 

G.S.  lost 

D. T.  Sherman 
A.D.  Graham 
J.M.  Grandage 
A.S.  Machin 
M.  Stimson 

P.  Piperias 

Materials  Research  Laboratory 
Director/Library 

Defence  Science  &  Technology  Organisation  -  Salisbury 
Library 

Navy  Office 

Navy  Scientific  Adviser  (3  Copies  Doc  Data  Sheet  Only) 

Director  of  Aircraft  Engineering  -  Navy 

Army  Office 

Scientific  Adviser  -  Army  (Doc  Data  Sheet  Only) 


Air  Force  Office 

Air  Force  Scientific  Adviser 
Engineering  Branch  Library 
Director  General  Engineering  -  Air  Force 
HQ  Logistics  Command  (DGELS) 

Statutory  and  State  Authorities  and  Industry 

Aero-Space  Technologies  Australia,  Systems  Division  Librarian 

Ansett  Airlines  of  Australia,  Library 

Australian  Airlines,  Library 

Qantas  Airways  Limited 

Gvil  Aviation  Authority 

Hawker  de  Havilland  Aust  Pty  Ltd,  Viaoria,  Library 
Hawker  de  Havilland  Aust  Pty  Ltd,  Bankstown,  Library 

Universities  and  Colleges 
Adelaide 

Professor,  Mechanical  Engineering 

Melbourne 

Engineering  Library 

Monash 

Head,  Materials  Engineering 

Newcastle 

Library 

Professor  R.  Tclfer,  Institute  of  Aviation 
Sydney 

Engineering  Library 
NSW 

Head,  Mechanical  Engineering 

Library,  Australian  Defence  Force  Academy 

Tasmania 

Engineering  Library 

Western  Australia 
Library 

Head,  Mechanical  Engineering 

RMIT 

Library 

Mr  M.L.  Scott,  Aerospace  Engineering 


CANADA 

CAARC  Coordinator  Structures 
NRC 

J.H.  Parkin  Branch  (Aeronautical  &  Mechanical  Engineering  Library) 

Universities  and  Colleges 
Toronto 

Institute  for  Aerospace  Studies 

FRANCE 

ONERA 

INDIA 

CAARC  Coordinator  Structures 

National  Aeronautical  Laboratory,  Infomnation  Centre 

INTERNATIONAL  COMMITTEE  ON  AERONAUTICAL  FATIGUE 
per  Australian  ICAF  Representative  (27  copies) 

ISRAEL 

Technion-Israel  Institute  of  Technology 
JAPAN 

National  Research  Institute  for  Metals,  Fatigue  Testing  Div. 
NETHERLANDS 

National  Aerospace  Laboratory  (NLR),  Library 
NEW  ZEALAND 

Defence  Scientific  Establishment,  Library 

Universities 

Canterbury 

Head,  Chemical  Engineering 


SINGAPORE 

Director,  Defence  Materials  Organisation 
SWEDEN 

Aeronautical  Research  Institute,  Library 
SWITZERLAND 

F+W  (Swiss  Federal  Aircraft  Factory) 


UNITED  KINGDOM 

Defence  Research  Agency  (Aerospace) 

Famborough,  Head  Materials  Department 
British  Library,  Document  Supply  Centre 
CAARC  Coordinator,  Structures 
British  Aerospace 

Kingston-upon-Thames,  Library 
Hatfield-Chester  Division,  Library 

Universitigs  and-Collggei 

London 

Head,  Aero  Engineering 
Manchester 

Professor,  Applied  Mathematics 
Head,  Dept  of  Engineering  (Aeronautical) 

Cranfield  Inst,  of  Technology 
Library 

Imperial  College 

Aeronautics  Library 

UNITED  STATES  OF  AMERICA 

NASA  Scientific  and  Technical  Information  Facility 
Boeing  Commercial  Airplanes 
Lockheed  California  Company 
Lockheed  Georgia 
McDonnell  Aircraft  Company 
Library 

Mr  J.  Young,  FI 8  Engineering  Dept. 

Northrop  Corporation,  Mr  J.  Knight,  FI  8  Structural  Analysis  Org. 

Universities  and  Colleges 

Florida 

Aero  Engineering  Department 
Head,  Engineering  Sciences 


SPARES  (10  COPIES) 


TOTAL  (126  COPIES) 


AL  140 


DCPARIMENT  OF  DEFACE 


< 


DOCUMENT  CONTROL  DATA 


PACE  CLASSmCATlON 

UNCLASSIFIED 


PRIVACY  MARKINC 


la.  AR  NUMBER 

AR-006-626 

lb.  ESTABUSHMENTNUMBER 

ARL-STRUC-R-i44 

2.  DOCUMENT  DATE 

SEPTEMBER  199 

1 

3.  TASK  NUMBER 

AIR  89/079 

A.  TITLE 

CRITICAL  APPRAISAL  OF  THE 

McDonnell  douglas  closure 

MODEL  FOR  PREDICTING  FATIGUE 
CRACK  GROWTH 

5.  ffiCURTTY  CLASSmCATION 

[PLACE  APTODPRIATE  CLASSZFlCA-nON 
tNBQ3qS]  IE.  SECRET  (SL  CONF.  (O 
RESTRICTED  [RJ.  UNCLASSIFIED  (U)  ]. 

6.  NO.  PACES 

85 

7.  NO.  REFS. 

19 

u  u 

U 

DOCUMENT  'TTILE  ABSTRACT 

8.  AUTHORIS) 

MJ.  BOS 

0.  DOWNCRADING/DELIMmNC  INSTRUCnONS 

Not  applicable 

10.  CORPORATE  AUTHOR  AND  ADDRESS 

AERONAUTICAL  RESEARCH  LABORATORY 

506  LORIMER  STREET 

FISHERMENS  BEND  VIC  3207 

n.  OFFTCE/POSmON  RESPONSIBLE  FOR; 

SI«NSOR  ^^AAF 

SECURITY 

DOWNGRADINC 

APPROVAL  DARL 

12.  SECONDARY  DISTRIBUTION  (OF  THIS  DOCUMENT 


Approved  for  public  release 

OVERSEAS  ENQUmiES  OUTSIDE  STATED  LIMITATIONS  SHOULD  BE  REFERRED  THROUGH  DSTIC,  ADMINISTRATIVE  SERVICES  BRANCH, 
DEPARTMENT  OF  DEFENCE.  AN2AC  PARK  WEST  OFFICES.  ACT  260 1 


13«.  THIS  DOCUMENT  MAY  BE  ANNOUNCED  IN  CATALOGUES  AND  AWARENESS  SERVICES  AVAU-ABLE  TO _ 

No  limitations 


I3b.  CITATION  FOR  OTHER  PURPOSES  (IE.  CASUAL 
ANNOUNCEMENP  MAY  BE 


UNRESTRICTED  OR 


□  AS  FOR  13i. 


14.  DESCRIPTORS 

IS.  DISCAT  SUBJECT 

Crack  propagation  modeling 

CATEGORIES 

Flight  loads 

0103 

Aircraft  structures 

110601 

ABSTRACr 

;  The  McDonnell  Douglas  fatigue  crack  growth  model  CNTKM8  has  been  used  to  predict  crack  growth  in 
\  CCT-sprc.-nens  with  though  cracks  under  four  different  load  sequences,  including  a  flight-by-flight 
sequence,  ana  f^o  stress  levels.  The  predictions  were  compared  against  experimental  data.  The  results 
ranged  from  fairly  good  for  the  more-structured  .sequences,  in  rather  poor  for  the  more-random  sequences. 
Furthermore,  the  experime.ntai  t'cru!  of  the  structured  .seourncr.s  giving  longer  crack  growth  lives  than  the 
random  sequences  was  not  predicted  correctly.  An  enhanced  modrl  was  proposed,  which  waj  also  tested. 
The  predictions  of  this  model  were  mere  consistent  with  the  experimental  data,  and  the  above-mentioned 
experimental  trend  was  also  predicted  correctly. 

This  report  contains  a  detailed  description  of  both  models,  including  the  cquadons  that  are  relevant  for  the 
case  of  a  CCT-specimen  with  a  through  crack,  and  of  the  tests  that  they  were  subjected  to. 


nwE  ciAasmcATioN 

UNCl>ASSIFIED 


PRIVACY  MAMONO 

ms  PAOS  IS  TO  BE  USED  TO  RECORD  INFORMATION  VWOCH  IS  REQUINED  BY  THE  EOTABUSHMENT  FOR  flB  OWN  USE  BUT  WinCH  WILL  NOT 
SC  ADOEDTO-niE  OIST1S  DATA  UNLESS  STECinCAaY  REQUESIEO. 


IS  ABSTRACT  ICONT1. 


IT.  IMPRINT 


AERONAUTICAL  RESEARCH  LABORATORY*  MELBOURNE 


1  18.  DOCUMENT  SCRlEfliUmNiafBCR 

10.  COOT  CODE 

20.  TYPE  OF  REPORTTAND  PERIOD  COVERED 

Aircraft  Structures  Report  444 

25205A 

7\.  COMJVI^ f>nOCRAM5 USED 


McAir  crack  closure  CNTKM8 
ARL  crack  closure  ARLECC 


;  73  rjCAinjsiMyrniKmrai 


5  1 


